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1.0  INTRODUCTION 


In  either  a  deterministic  or  probabilistic  seismic  hazard  evaluation  an 
essential  element  is  a  description  of  the  variability  of  ground  motion 
parameters  with  distance  (or  depth)  and  earthquake  size  or  magnitude. 

Once  a  distance  and  design  basis  earthquake  is  decided  upon,  the  region 
specific  attenuation  relation  is  then  utilized  as  an  estimator  to  predict, 
with  associated  variance,  the  ground  motion  parameters  to  be  expected  at 
the  site.  If  attenuation  relations  are  available  which  parameterize  peak 


values  and  frequency  content  through  the  response  spectra,  then  the 
exposure  is  reasonably  well  specified.  If  however,  only  peak  values  are 
constrained  by  the  attenuation  relation,  these  may  be  utilized  to  scale  an 
assumed  frequency  dependence  via  an  adopted  response  spectrum  shape. 

While  not  as  satisfying  as  a  region  specific  response  spectrum  shape,  the 
adopted  shapes  are  based  upon  many  observations  and  can  be  adjusted  to 
relect  the  magnitude  range  contributing  to  the  seismic  hazards. 

The  representation  of  the  design  ground  motions  by  a  response  spectrum  is 
sufficient  for  the  seismic  design  evaluations  of  most  engineered 
structures  (e.g.  commercial  buildings,  hospitals,  etc.).  However,  for 
large,  complex  and/or  critical  facilities,  time  history  analyses  are  often 
necessary,  especially  when  there  is  significant  non-linear  structural 
response  (e.g.  earth  dams,  offshore  platforms,  etc.).  For  these 
evaluations,  one  or  more  accelerograms  are  selected  whose  response 
spectrum  match  the  design  spectrum  in  some  average  sense. 

In  regions  of  high  seismicity  rates  which  have  had  established  strong 
motion  instrumentation  programs  for  a  period  of  time,  a  sufficient  data 
base  of  observations  may  be  available  to  provide  representati ve 
accelerograms  and  to  constrain  regression  analyses  for  peak  values  and 
response  spectra.  While  the  range  of  magnitudes  and  distances  is  always 
less  than  ideal,  definite  trends  may  be  characterized  statistically  so 


that  extrapolations  may  be  made  with  some  level  of  confidence  in  not  only 
the  values  themselves  but  also  the  degree  of  conservatism  as  well. 


Time  histories,  which  are  consistent  in  amplitude,  duration,  and  frequency 
content  with  observational  data  are  generated  through  filtering  noise 
samples  (Nau  et  al.,  1982)  or  by  splicing  together  selected  portions  of 
observed  time  histories.  Also,  response  spectral  matching  techniques  may 
be  applied  to  empirical  data  which  scale,  through  filters,  the  time 
history  (Tsai,  1972).  This  results  in  a  time  history  whose  response 
spectra  is  close  to  the  design  response  spectrum. 

Although  these  approaches  retain  the  frequency,  amplitude,  and  duration 
characteristics  in  a  statistical  sense,  they  usually  produce  time 
histories  which  could  not  have  resulted  from  any  earth  process.  The 
resulting  velocity  and  displacement  time  histories  from  the  artificially 
generated  accel erograms  are  often  inconsistent  with  observed  motions,  a 
fact  that  may  be  essential  in  response  calculations  for  longer  period 
structures.  These  approaches  have  evolved  through  necessity  and  employ 
little  knowledge  of  earthquake  source  processes  or  wave  propagation 
physics. 

At  the  other  extreme  in  ground  motion  prediction,  detailed  descriptions  of 
source  properties  (e.g.  location,  direction,  rupture  velocity, 
distribution  of  asperities  or  barriers,  and  rise  times)  and  propagation 
path  parameters  are  utilized  to  generate  synthetic  time  histories  (Heaton 
and  Helmberger,  1978,  Apsel  et  al.,  1983).  These  models  are  essential  to 
understand  past  seismic  events  and  they  can  be  utilized  to  predict  future 
ground  motions.  However,  they  require  the  specifications  of  source 
details  and  path  effects  which  can  have  profound  effects  on  the  predicted 
ground  motion.  Generally,  a  range  in  parameters  is  utilized  in  a 
sensitivity  analysis  which  translates  into  uncertainties  in  the  predicted 
motion. 

The  modeling  approach  has  recently  attempted  to  accommodate,  in  a  natural 
way,  some  of  the  stochastic  aspects  of  source  processes  and  path  effects 
by  utilizing  small,  well  recorded  events  as  a  basis  to  construct  time 
histories  (Hartzell,  1978).  The  observed  events  are  scaled  and  summed, 
utilizing  source  physics,  to  model  time  histories  for  large  events  (Hadley 


and  Helmberger,  1980;  Kanamori,  1979).  This  hybrid  approach  is  very 
attractive  in  that  aspects  of  both  source  and  path  complexities  are 
incorporated.  This  reduces  somewhat  the  arbitrary  nature  in  parameter 
specification  and  allows  more  definitive  calibrations  with  large  events  to 
help  constrain  the  remaining  variables. 

An  extremely  important  application  of  the  hybrid  approach  is  its  utility 
in  providing  guidelines  in  extrapolations.  Estimates  of  ground  motion 
parameters  using  various  empirical  models  (Donovan,  1973;  Idriss,  1978; 
Joyner  and  Boore,  1981;  Campbell,  1981;  Joyner  and  Boore,  1982;  Joyner  and 
Fumal ,  1984)  differ  little  in  regions  of  distance  and  magnitude  where  data 
are  abundant.  However,  in  regions  where  extrapolations  are  required,  close 
to  moderate  and  large  earthquakes,  the  differences  can  be  large.  The 
differences  are  due  primarily  to  the  mathematical  form  chosen  to 
parameterize  the  data.  In  these  regions,  the  careful  use  of  modeling 
which  is  based  upon  observational  data  can  be  a  powerful  tool  in  assessing 
the  conservatism  of  extrapolated  empirical  relations  (Hadley  and 
Helmberger,  1982).  The  modeling  is  calibrated  in  regions  where  data  exist 
by  parameter  variation  and  sensitivities  are  assessed.  A  level  of 
confidence  is  achieved,  at  least  by  the  modelers,  and  estimates  are  made 
on  ground  motion  which  may  guide  or  give  confidence  in  extrapolated 
empirical  relations. 

While  the  above  techniques  are  the  basic  tools  in  ground  motion 
estimation,  they  are  both  fundamentally  based  upon  observational  data. 

The  empirical  approach  obviously  cannot  be  utilized  with  few  strong  motion 
data,  while  the  analytical  and  hybrid  techniques  utilize  data  either 
directly  or  at  least  for  calibration  purposes  to  constrain  free 
parameters. 

In  many  areas  of  the  world,  however,  there  exists  active  tectonism  and 
therefore  considerable  seismic  exposure  but  the  strong  motion  data  base  is 
extremely  sparse.  In  these  cases,  an  empirical  relation  based  upon  non¬ 
region  specific  data  may  simply  be  adopted.  If  time  histories  are 
required,  conventional  scaling  techniques  are  employed  generally  using  WUS 
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acceleration  data.  The  justification  primarily  relies  upon  arguments 
considering  similar  tectonic  environments,  crustal  structure,  style  of 
faulting,  depth  of  events,  observed  attentuation  of  shaking  intensity  and 
other  less  tangible  aspects  generally  referred  to  as  engineering 
judgement.  While  this  methodology  is  an  accepted  practice  by  necessity, 
an  approach  is  needed  which  does  not  rely  upon  strong  motion  data  for 
calibration  but  rather  only  for  confirmation. 
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2.0  APPROACH 


Another  approach,  which  utilizes  simple  earthquake  source  theory  and  wave 
propagation  physics,  has  recently  demonstrated  great  promise  in 
circumstances  such  as  these  (Boore,  1983;  Atkinson,  1984,  McGuire  et  al., 
1984).  The  basic  advantage  lies  in  using  weak  motion  indirectly  to 
predict  strong  motion.  The  weak  motion  is  due  to  small  magnitude 
(generally  less  than  Mw  =  5)  local  or  regional  events.  These  data,  from 
both  analogue  and  digital  recordings,  are  used  to  estimate  region  specifi 
source  and  wave  propagation  parameters.  Those  parameters  are  then  input 
to  the  source  and  propagation  models  which  predict  the  motion  due  to 
events  at  distances  and  from  source  sizes  for  which  no  data  exist.  This 
approach  employs  random  vibration  theory  (RVT)  applied  to  a  Brune  source 
spectrum  (Brune;  1970,  1971)  to  characterize  strong  ground  motion.  While 
this  approach  is  not  perfect,  and  some  problems  remain,  we  believe  that, 
based  upon  its  recent  success  (Boore,  1983;  Atkinson,  1984,  McGuire  et 
al.,  1984)  and  our  preliminary  results,  it  will  prove  to  be  quite  useful. 

2.1  Model 

The  formalism  employed  to  develop  an  attenuation  relation  with  little  or 
no  strong  motion  data  is  based  upon  a  simple  theoretical  model  of  the 
earthquake  process  and  wave  propagation  physics.  The  theoretical  basis 
was  developed  and  calibrated  with  observed  data  by  Hanks  and  McGuire 
(Hanks,  1979;  McGuire  and  Hanks,  1980;  Hanks  and  McGuire,  1981).  They 
used  the  Brune  (1970,  1971)  spectrum  to  model  root  mean  square  (RMS) 
acceleration  as  a  function  of  magnitude  and  distance  for  stiff  (rock) 
sites.  To  model  peak  values,  they  utilized  results  from  random  vibration 
theory  to  relate  the  RMS  predictions  to  maximum  values.  Boore  (1983)  has 
extended  the  range  of  applications  to  include  predictions  of  peak 
horizontal  velocities,  Wood-Anderson  seismographic  response,  and  response 
spectra.  Boore  has  also  used  this  approach  to  generate  synthetic 
acceleration  time  histories  by  employing  random  sequences  whose  spectra 
match  the  predicted  Brune  spectrum.  The  method  leads  to  results  which 
reproduce  the  empirical  dependences  of  peak  acceleration,  peak  velocity. 
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reproduce  the  empirical  dependences  of  peak  acceleration,  peak  velocity, 
and  pseudovelocity  response  spectral  amplitudes  on  moment  magnitude  at 
close  distances.  McGuire  et  al.  (1984)  further  confirmed  the  results  of 
the  RVT  technique  by  demonstrating  the  close  agreement  between  observed 
and  predicted  response  spectral  values.  They  also  showed  good  agreement 
between  the  calculated  Brune  spectrum  and  the  acceleration  spectral 
density  of  empirical  data  for  several  discrete  frequencies.  In  order  to 
extend  this  approach  to  large  distances,  the  model  was  modified  to 
incorporate  surface  wave  effects  (Atkinson,  1984).  At  distances  greater 
than  about  two  crustal  thickness  (Herrmann,  1985),  surface  waves  rather 
than  shear  waves,  will  carry  the  peak  ground  motion. 


The  technique  developed  in  this  study  employs  aspects  of  the  above 
developments  to  predict  peak  values  (acceleration  and  velocity)  and 
response  spectra.  However  it  differs  significantly  from  other  techniques 
in  the  manner  ir  which  time  histories  are  generated.  These  time  histories 
may  be  regarded  as  semi -empi rical  in  that  they  use  the  Brune  spectrum  as  a 
modulus  but  an  observed  phase  to  generate  the  complex  spectrum.  The 
advantage  of  this  technique  is  that  the  non-stationarity,  randomness,  and 
change  in  frequency  with  time  is  incorporated  in  a  natural  way. 
Integrations  to  velocity  and  displacement  are  then  also  more  realistic. 
The  basic  assumption  with  this  technique  is  that  the  region-specific 
source  and  wave  propagation  parameters  are  reflected  primarily  in  the 
spectral  modulus.  The  phase  spectrum  accounts  for  the  multipath  effects 
and  surface  wave  contributions. 


Since  the  RVT  technique,  employing  a  Brune  spectrum,  has  demonstrated 
success  in  modeling  both  Western  and  Eastern  United  States  data,  the  RVT 
peak  values  (acceleration  and  velocity)  are  used  to  scale  the  time  domain 
simulations.  This  is  necessary  since  we  are  generally  combining  a  modulus 
and  a  phase  which  are  somewhat  incompatible.  This  arises  because  the 
phase  is  from  a  time  history  with  mixed  phase  properties  while  the  Brune 
spectrum  is  smooth  (Pilant  and  Knopoff,  1970). 
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In  order  to  preserve  the  magnitude  dependency  of  strong  motion  duration 
characteristics  (Oobry  et  al.,  1978)  it  is  necessary  to  use  a  phase  from  a 
record  of  approximately  the  same  magnitude  and  distance  as  the  design 
event.  This  arises  because  the  phase  spectrum  determines  how  the  energy 
is  distributed  in  time.  The  magnitude  similarity  is  about  _+  one  half  unit 
of  magnitude.  The  distance  requirement  on  the  phase  is  less 
restrictive.  Generally,  we  have  found  phases  extracted  from  records 
within  20  to  25  km  can  be  used  for  simulations  from  10  to  30  km.  Phases 
from  records  at  40  to  50  km  are  used  for  distances  greater  than  30  km. 

In  cases  where  a  time  history  is  required  when  a  target  response  spectrum 
is  specified,  an  initial  Brune  spectrum  is  generated  based  upon  magnitude, 
stress  drop,  and  distance  in  addition  to  the  region  specific  wave 
propagation  parameters.  The  Brune  spectrum  is  then  scaled  by  taking 
ratios  of  either,  the  RVT  response  spectra  or  a  time  domain  response 
spectrum  calculation,  to  the  target  response  spectrum.  The  time  domain 
response  spectrum  is  calculated  from  the  synthesized  time  history.  The 
peak  acceleration  of  the  resultant  time  history  is  generally  then  scaled 
to  the  design  value  (which  must  be  consistent  with  the  target  response 
Spectrum) . 


3.0  OBJECTIVE 


The  original  objective  of  this  work  was  to  develop  a  technique  to  generate 
a  realistic  time  history  whose  response  spectrum  is  compatible  with  a 
specified  response  spectrum.  The  plan  was  to  use  an  observed  time  history 
and  scale  the  modulus  of  its  Fourier  spectral  density  by  the  ratio  of  the 
observed  response  spectra  to  the  target.  The  original  phase  would  then  be 
added  to  produce  the  desired  time  history.  The  process  would  be  iterated 
until  the  response  spectra  of  the  scaled  time  history  was  close  to  the 
target  response  spectrum.  However,  in  developing  the  code  it  was  decided 
to  abandon  the  observed  modulus  since  the  RVT  technique  had  recently 
demonstrated  such  good  results  by  employing  the  simple  Brune  spectrum. 

The  most  convincing  evidence  came  from  McGuire  et  al.,  (1984).  They 
showed  the  effectiveness  of  the  Brune  spectra  in  predicting  Fourier 
spectral  density  and  response  spectral  values  at  close  distances  for  several 
frequencies.  By  using  a  combined  approach,  that  is  a  Brune  modulus  with  an 
observed  phase,  it  is  possible  to  also  incorporate  all  the  predictive  power 
of  the  RVT  methodology.  That  is,  in  one  code,  the  ability  to  predict  peak 
acceleration,  peak  velocity,  and  response  spectra  in  addition  to  generate 
realistic  acceleration,  velocity,  and  displacement  time  histories  has  been 
incorporated.  This  can  be  done  for  magnitudes  from  about  4  to  7  1/2  and 
from  distances  from  10  km  to  over  100  km.  In  addition  a  target  response 
spectrum  may  be  input  to  the  code  and  an  acceleration  time  history  can 
then  be  generated  whose  response  spectrum  closely  matches  the  target 
response  spectrum. 
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4.0  RESULTS 


Since  the  code,  as  it  is  presently  configured,  operates  in  two  basic  modes 
depending  upon  whether  or  not  a  target  response  spectrum  has  been 
specified,  the  results  section  will  be  divided  accordingly.  The  first 
section  will  address  the  predictive  capability  in  terms  of  peak  values  for 
WUS  and  EUS  tectonic  environments.  The  time  history  synthesis  will  also 
be  demonstrated.  The  second  section  will  present  the  scaling  results  when 
an  input  target  response  spectrum  is  specified. 

4.1  Prediction  and  Synthesis  Capability 


In  order  to  demonstrate  the  predictive  characteristics  of  the  code, 
results  for  both  magnitude  and  distance  scaling  will  be  presented. 
Magnitude  scaling  will  include  comparisons  of  peak  acceleration  and  peak 
velocity  at  close  distances  with  empirical  relations  for  WUS  and  EUS 
tectonic  environments.  Distance  scaling  will  involve  predicted  peak 
acceleration  vs  distance  for  a  single  magnitude  compared  with  empirical 
data.  Source  and  wave  propagation  parameters  used  for  WUS  and  EUS 
predictions  are  shown  in  Table  1. 


?w\ 


Figure  1  demonstrates  the  magnitude  scaling  for  close  distances  (R<15  km) 
compared  to  three  other  attenuation  relations:  Joyner-Boore  (1982)  at  50% 
level,  Seed  and  Schnable  (1980),  and  Donovan  (1973).  The  bars  at 
magnitudes  5  and  7  1/2  represent  the  scatter  in  the  data  base  used  by 
Joyner  and  Boore  (1981).  The  lower  magnitude  variances  are  for  the 
Oroville  aftershocks  combining  rock  and  soil  sites  (Shakal  and  Berneuter, 
1980).  All  data  are  for  distances  of  a  few  km  to  15  km.  Both  the  Joyner- 
Boore  (1982)  and  Seed  and  Schnable  (1980)  curves  are  plotted  for  magnitude 
ranges  over  which  the  authors  indicated  the  relations  are  valid.  These 
relations  are  based  upon  WUS  data.  The  Donovan  (1973)  curve,  which  is 
based  upon  world  wide  data,  is  plotted  for  the  entire  range  since  no 
discussion  is  given  outlining  the  range  of  its  validity.  It  should  also 
be  pointed  out  that  there  are  differences  in  definitions  of  distance  which 
are  significant  at  these  close  ranges.  The  calculated  relation  is  based 
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upon  a  hypocentral  distance  of  10  km  while  the  Joyner-Boore  (1982)  curve 
uses  the  closet  distance  to  the  surface  projection  of  the  fault.  Seed- 
Schnable  (1980)  and  Donovan  (1973)  employ  a  "closest"  distance.  These 
differences  in  definition  can  easily  result  in  differences  of  0.2  to  0.3 
log  units  at  these  close  distances  (Hanks  and  McGuire,  1981).  The  RVT 
predictions,  which  include  the  near  site  amplification  factors  (Boore, 
1985),  appears  somewhat  high  but  is  generally  within  the  85%  Joyner  Boore 
(1982)  curve  (0.27  log  units  higher  than  the  curve  shown  in  figure  1). 
Also,  for  the  higher  magnitudes,  10  km  is  well  within  the  source  dimension 
which  violates  the  far-field  assumption  in  the  Brune  (1970,  1971) 
theory.  With  these  considerations,  we  are  favorably  impressed  with  the 
resul ts . 

Peak  particle  velocities  vs  moment  magnitude  are  plotted  in  Figure  2  for 
WUS  parameters.  Also  shown  is  the  empirical  curve  of  Joyner  and  Boore 
(1982).  Here  the  agreement  between  the  calculated  values  and  the 
empirical  relation  is  quite  good  over  the  entire  range  in  magnitudes. 

Distance  scaling  of  peak  acceleration  is  shown  in  Figure  3  for  a  moment 
magnitude  of  7  at  a  focal  depth  of  10  km.  The  data  are  from  the  1979 
Imperial  Valley  event  (M$  =  6.8;  Idriss,  1  983).  The  median  and  _+  1  (X 
curves  are  from  a  fit  to  the  data  (Idriss,  1983).  The  calculated  values 
(open  symbols)  generally  fall  within  the  data  indicating  a  realistic 
distance  seal ing. 

Eastern  United  States  scaling  is  shown  in  Figure  4.  In  this  case,  the 
near  site  amplification  factors  are  not  employed  since  the  velocity 
gradients  in  the  top  1  km  or  so  are  much  less  in  the  EUS.  The  plot 
compares  predicted  EUS  peak  acceleration  values  vs  magnitude  (mb)  to 
several  other  relations.  The  figure  was  taken  from  Atkinson  (1984).  The 
predicted  values  (open  symbols)  are  generally  within  the  range  of  the 
other  relations.  This  indicates  that  the  model  agrees  with  other 
predictions  of  magnitude  scaling  at  close  distances  for  the  EUS  tectonic 
environment.  Figure  5  shows  the  same  results  for  peak  velocity.  Here  a 
much  wider  range  is  shown  by  the  relations  but  again  the  model  predictions 
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are  generally  consistent  with  those  obtained  using  other  relations. 

Distance  scaling  for  the  peak  acceleration  in  the  EUS  is  shown  in  Figure 
6.  The  data  are  from  the  Miramichi,  St.  Lawrence,  New  Hampshire  areas  and 
are  scaled  to  mb  =  5.0  (see  figure).  The  line  is  from  Atkinson  (1984)  and 
represents  her  RVT  results  for  Eastern  Canada.  The  model  predictions 
(open  symbols)  are  for  a  source  depth  of  4  km  and  lie  well  within  the 
scatter  of  the  data.  The  figure  was  taken  from  Atkinson  (1984). 

An  important  consideration  in  the  analysis  of  degrading  systems  such  as 
liquifiable  soil  deposits  and  saturated  earth  embankments  is  the  number  of 
cycles  of  loading.  This  is  related  to  the  time  history  through  measures 
of  duration.  In  order  to  assure  that  the  synthesized  time  histories 
reflect  an  appropriate  increased  duration  with  magnitude,  the  significant 
durations  were  monitored  for  moment  magnitude  4  to  7  and  at  hypocentral 
ranges  of  10  and  50  km.  The  results  are  shown  in  Figure  7  along  with  the 
empirical  curve  of  Dobry  et  al.  (1978)  for  the  5%  to  95%,  Arias  intensity 
at  close  distances.  The  Dobry  et  al.  curve  is  shown  over  the  range  of 
their  data  and  the  synthetic  results  agree  well  for  magnitudes  6  to  7. 

The  agreements  degrades  for  lower  magnitudes  but  the  scatter  in  this  type 
of  data  is  large.  In  general,  the  synthetic  results  shown  the  expected 
trend  with  magnitude  and  distance. 

The  corner  period  is  also  shown  vs  moment  magnitude  in  Figure  7.  This 
relation  (Table  1)  is  nearly  within  _+  2  and  demonstrates  that  the 
inverse  corner  frequency  (source  duration)  can  be  a  good  measure  of  the 
lower  bound  estimates  of  strong  ground  motion  duration  at  close  distances. 

The  time  domain  synthesis  is  shown  in  Figure  Set  8  for  both  Western  and 
Eastern  United  States  scaling  (Table  1).  Acceleration,  velocity,  and 
displacement  time  histories  are  shown  for  a  magnitude  7  (Mw  for  WUS  and  mb 
for  EUS)  event  at  10  km.  Peak  time  domain  accelerations  and  velocities 
are  scaled  to  the  corresponding  RVT  peaks  while  the  displacements  are 
integrated  from  the  scaled  accelerations.  For  many  applications  is  site 
response  analysis,  it  is  desirable  to  produce  time  histories  and  response 


spectra  at  shallow  depths.  This  would  correspond  to  embedded  structures 
with  foundations  within  a  soft  or  stiff  site.  In  order  to  accommodate 
this  feature,  the  code  includes  a  single  layer  site  response  transfer 
function  (Section  6).  While  most  sites  have  material  properties  which 
vary  with  depth,  the  overall  effects  of  burial  may  be  represented  by  a 
single  layer  with  correctly  averaged  properties.  The  results  for  a 
receiver  at  a  depth  of  20  m  within  a  40  m  thick  soil  sites  (Table  3)  are 
shown  in  Figure  Set  8  for  WUS  parameters.  The  respective  response  spectra 
are  also  shown  in  Figure  Set  8.  The  depth  dependent  spectral  node  at 
around  4  Hz  is  clearly  visisble  in  the  at-depth  response  spectrum. 

An  example  of  the  RVT  response  spectra  is  shown  in  Figure  Set  9.  The 
upper  solid  line  in  plot  (A)  is  the  pseudo  relative  velocity  response 
spectrum  for  a  moment  magnitude  6  event  at  a  distance  of  10  km  calculated 
using  the  RVT  formulation.  The  scaling  is  for  the  WUS  (Table  1)  and  the 
open  circles  were  taken  from  Boore  (1983)  and  represent  regression  fits  to 
WUS  data.  The  agreement  is  reasonably  good  and  well  within  the 
variability  of  the  data.  We  should  also  recall  that  the  RVT  calculations 
represent  expected  values  (Udwadia  and  Trifunac,  1974).  The  5  to  95% 
confidence  limits  depend  upon  N,  number  of  peaks)  but  in  general  the  range 
is  about  50%  above  and  below  the  expected  value  (Udwadia  and  Trifunac, 
1974).  The  kinks  at  the  long  period  end  of  the  RVT  response  spectrum  are 
due  to  N  changing  with  frequency.  At  long  periods  N  is  small  so  a  change 
in  one  unit  has  a  substantial  effect  on  the  peak  to  RMS  ratio. 

4.2  Response  Spectral  Scaling  Capability 

In  order  to  demonstrate  the  response  spectral  scaling  capability  for  the 
code,  a  design  spectrum  from  a  previous  Woodward-Clyde  project  has  been 
utilized  in  order  to  compare  results  of  the  present  formulation  with 
current  practice. 

The  response  spectral  scaling  approach  taken  in  this  code  employs  both  RVT 
and  time  domain  response  spectra  calculations.  Since  RVT  response  spectra 
calculations  are  less  costly  than  time  domain  evaluations,  the  first 


iterations  (1  or  2)  are  done  with  this  technique.  This  also  provides  for 
a  more  stable  convergence  since,  during  the  first  few  iterations,  the 
Brune  spectra  are  perturbed  the  greatest,  therefore  using  extremely  smooth 
RVT  response  spectra  results  in  non-oscil 1 atory  scaling  factors.  Since 
the  RVT  response  spectra  are  expected  values,  the  response  spectra  calcu¬ 
lated  from  the  synthesized  time  history  may  depart  from  these  estimates  by 
as  much  as  ±  50%.  This  variability  is  due  to  the  use  of  an  observed  phase 
spectrum  and  reflects  path  and  site  effects.  To  correct  for  this,  the 
final  iterations  (1  or  2)  are  done  employing  response  spectrums  calculated 
from  the  synthesized  time  histories.  Direct  integration  of  the  oscillator 
equation  is  performed  (Nigam  and  Jennings,  1968).  For  oscillator  periods 
less  than  10  times  the  sample  interval,  the  acceleration  time  histories 
are  linearly  interpolated.  This  was  done  to  provide  more  accurate  high 
frequency  response  calculations. 

Figure  Set  9  shows  the  target  or  design  response  spectrum  and  the  scaling 
results.  The  design  event  is  taken  to  be  a  WUS  earthquake  with  a  moment 
magnitude  of  6  at  a  distance  of  20  km.  The  peak  acceleration  is  specified 
at  0.125  g.  The  high  frequency  limit  (50  Hz)  of  the  target  PRSV  gives  a 
peak  acceleration  of  0.123  g  so  the  design  peak  acceleration  and  the  high 
frequency  limit  of  the  PSRV  are  consistent.  Plot  (A)  shows  the  design 
PSRV  and  the  initial  RVT  response  spectrum.  Also  shown  are  the  results 
after  the  first  two  iterations  employing  RVT  spectrum  calculations  and  the 
final  two  iterations  using  the  time  domain  response  spectra 
calculations.  The  convergence  is  remarkably  fast  and  the  fit  is  good. 

Plot  (B)  shows  the  Fourier  spectral  density  after  each  set  of  two 
iterations.  The  oscillations  introduced  into  the  Fourier  spectra  when 
using  time  domain  response  spectra  calculations  for  scaling  are  readily 
apparent.  It  is  interesting  to  note  that  these  oscillations  are  due  to 
the  randomness  of  the  observed  phase  since  the  RVT  scaled  Brune  spectra 
are  quite  smooth  (dotted  line  in  Plot  (B)).  Thus  the  randomness  in  the 
observed  phase  is  indirectly  coupled  into  modulus  through  the  scaling 
factors.  The  Fourier  spectra  then  takes  on  a  more  realistic  appearance. 

In  order  to  provide  a  meaningful  integration  to  velocity  and  displacement, 
the  final  Fourier  spectra  are  band-pass  filtered  with  five-pole  causal 
Butterworth  filters.  In  this  case  the  high-pass  corner  was  0.25  Hz  and 


the  low-pass  corner  is  23  Hz.  The  results  of  the  filtering  on  the 
response  spectra  are  shown  in  plot  (C).  The  effects  of  the  high-pass 
filter  propagate  up  to  about  0.7  Hz. 

The  acceleration,  velocity,  and  displacement  time  histories  are  shown  in 
plots  (d)  of  Figure  Set  9.  The  time  histories  appear  very  realistic  with 
little  wrap-around  problems.  The  acceleration,  with  a  peak  of  0.116  g, 
displays  typical  non-stationarity  and  change  in  frequency  with  time.  The 
integrations  to  velocity  and  displacement  display  reasonable  values  and 
shapes. 

In  order  to  produce  the  design  peak  acceleration  of  0.125  g,  the 
acceleration  time  history  was  normalized  to  this  value  after  the  last 
iteration.  The  resulting  response  spectra  and  time  history  are  shown  in 
Figure  Set  10.  Plot  (a)  shows  the  response  spectra  which  is  now  slightly 
above  the  target  response  spectra.  This  arises  because  the  normalization 
process  to  scale  the  peak  acceleration  to  0.125  g  scales  the  entire  time 
history  resulting  in  a  perturbed  response  spectrum  throughout  the  entire 
bandwidth  of  oscillator  frequencies.  A  more  proper  manner  of 
accomplishing  this  would  be  to  perturb  the  Fourier  spectra  in  the 
frequency  range  which  controls  the  time  domain  peak.  From  Figure  Set  9, 
plot  (B),  this  can  be  seen  to  occur  around  2  Hz.  Nevertheless,  we 
consider  the  fit  quite  acceptable. 

Figure  Set  11  shows  the  results,  for  the  same  target  response  spectra  and 
peak  acceleration  employing  another  scaling  technique.  The  approach  used 
here  was  to  piece  together  elements  from  several  time  histories  until  the 
resulting  response  spectra  approximated  the  design  response  spectra.  At 
that  point,  spectrum  raising  and  lowering  techniques  were  applied  (Tsai, 
1972)  to  provide  the  match  shown  in  the  first  plot.  The  fit  is  quite  good 
but  the  time  histories  appear  somewhat  unrealistic.  The  velocities  and 
displacements  are  naturally  less  satisfying  than  the  acceleration  and  may 
represent  little  more  than  the  results  of  processing. 


5.0  CONCLUSION  AND  RECOMMENDATIONS 


5.1  Conclusion 

In  general,  we  are  favorably  impressed  with  the  results  of  the  code.  The 
use  of  RVT  methodology  in  the  prediction  aspect  for  peak  values  and 
response  spectra  appears  to  yield  appropriate  results  with  both  WUS  and 
EUS  magnitude  and  distance  scaling.  The  simple  and  straigthforward 
approach  to  time  domain  synthesis  produces  realistic  earthquake 
acceleration,  velocity,  and  displacement  records.  The  time  histories 
display  characteristic  non-stationarity  and  change  in  frequency  content 
with  time  without  resorting  to  envelope  functions  or  other  devices.  The 
acceleration  records  additionally  show  expected  increases  in  duration  with 
magnitude  and  distance. 

The  response  spectral  scaling  methodology  resulted  in  very  rapid  conver¬ 
gence  (2-4  iterations)  to  the  target  design  spectrum.  The  resulting 
spectrum  compatible  time  histories  again  displayed  very  realistic 
properties,  particularly  when  compared  to  results  employing  a  conventional 
technique.  Another  advantage  of  the  RASCAL  code  in  response  spectral 
scaling  applications  is  the  ease  of  use  and  cost  effectiveness.  Once  a 
design  spectrum  has  been  generated,  one  run  of  the  code  may  be  all  that  is 
required  versus  about  2  man  weeks  for  conventional  techniques. 


6.0  MATHEMATICAL  DEVELOPMENT 
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[  Following  McGuire  and  Hanks  (1980)  and  Boore  (1983)  the  modulus  of  the 

acceleration  spectral  density  for  a  Brune  (1970,  1971)  pulse  may  be 
written  as 

p  Wf-J  (?“>() 

rrtt)* 

where  0.78  accounts  for  the  free  surface  effect  (2),  participation  of 
energy  into  two  horizontal  components  (2-1/2),  and  an  average  radiation 
pattern  (0.5b). 

The  source  and  wave  propagation  parameters  are: 
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Tectonic  specific  scaling  is  achieved  through  the  above  region  dependent 
parameters.  Magnitude  dependence  may  be  incorporated  through  specifying  a 
moment-magnitude  relation.  For  WUS  the  moment  magnitude  is  generally 
used. 


log  M0  =  1.5  Mw  +  16.1 


(Hanks  and  Kanamori,  1979)  (5.4) 


while  for  these  simulations  the  following  EUS  mb  relation  was  used 


log  M0  =  1.75  mb  +  14.12 


( EPRI ,  1985)  (5.5) 


Random  vibration  theory  (RVT  is  implemented  by  specifying  the  relationship 
between  peak  values  (acceleration,  velocity,  oscillator  response  for 
response  spectra)  and  root  mean  square  (RMS)  values.  RMS  is  calculated  in 
the  frequency  domain  through  Parseval's  relation  (Aki  and  Richards,  1981) 
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where  T is  a  measure  of  the  duration  of  the  acceleration  time  history 
(Vanmarke  and  Lai,  1980). 

The  relationship  between  expected  value  of  the  peak  (E(,  ))  and  the  RMS 
value  (aRM^)  is  based  upon  the  work  of  Cartwright  and  Longuet-Higgens 
(1956)  and  Udwadia  and  Trifunac  (1974).  The  following  method  is  taken 
from  Boore  (1983,  1985)  who  applied  it  to  estimating  peak  acceleration, 
peak  velocity,  response  spectral  ordinates,  Wood-Anderson  response,  and  to 
evaluate  spectral  scaling  relations. 


The  following  development  is  in  terms  of  predicting  peak  acceleration  but 
the  application  to  peak  velocity  is  identical.  The  application  to 
prediction  of  response  spectra  is  also  demonstrated. 

The  expected  value  of  the  peak  acceleration  is  given  by 


AS 


where  a^g  is  calculated  through  Equation  (5.6)  utilizing  the  Brune 
spectra,  CN^  are  the  binominal  coefficients  (N!|lt  (N-l)l)  and  ?  is  a 
measure  of  the  bandwidth: 

T  '  ^  (5.8) 


The  mk  are  moments  of  the  energy  density  spectrum  given  by 
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The  limit  of  the  sum  (N)  is  defined  as  the  number  of  extrema  in  the 
acceleration  time  history  of  duration  *1*  which  is  taken  as  fc_1  (Hanks  and 
McGuire,  1981).  An  estimate  of  this  number  is  given  by 


N  -  2  Tt 


(5.  10) 
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where  f  is  the  predominant  frequency 


(5.11) 

The  minimum  value  of  N  is  set  at  2  (Boore,  1983)  while  if  N  is  greater 
than  20  an  asymptotic  expression  is  used  for  £(3^); 

F  *  at»t  £<?■*''*')  7  +  (7^7 ;  i  7 

(Clough  and  Penzine,  1975)  (5.12) 


where  N,  still  defined  by  Equation  (5.10),  is  the  number  of  zero  crossings 
in  the  time  *T*  .  In  this  case  the  predominant  frequency  (7)  is  defined  by 


(5.12) 


The  a^g  ig  given  by 
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from  Equations  (5.6)  and  (5.9) 


In  order  to  employ  the  RVT  technique  to  estimate  the  response  spectrum, 
the  oscillator  transfer  function  is  incorporated  into  Equation  (5.9).  The 
peak  value  from  Equation  (5.7)  or  Equation  (5.12)  is  then  the  response 
spectral  ordinate  for  a  particular  oscillator  damping  and  resonant 
frequency.  To  see  this  we  may  write  out  the  inhomogenous  oscillator 
Equation 
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(5.14) 


where  X  is  the  oscillator  displacement,  the  damping,  and  a(t)  the 
acceleration  time  history.  Taking  Fourier  transforms  and  multiplying  by 
fn2  <Cn  “  )  the  modulus  of  Equation  (5.14)  may  be  written  as 
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where  f  x  is  the  pseudo  acceleration  (Hudson,  1979).  The  psuedo 
acceleration  transfer  function  is  then 
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and  when  H  (fn,^  /  f)  is  included  in  the  intergrand  of  Equation  (5.9) , 
the  resultant  RVT  peak  is  the  pseudo  acceleration  spectral  ordinate  for  a 
given  >>  and  fn.  The  present  code  uses  Simpson's  rule  to  perform  the 
integrations . 

For  response  spectra  calculations,  a  modification  is  needed  to  the 
duration*"/"* •  This  arises  because  for  short  duration  time  histories,  the 
longer  period  oscillators  do  not  have  sufficient  time  to  build  up  their 
RMS  response.  Boore  and  Joyner  (1984)  have  an  empirical  correction  factor 
which  employs  an  equivalent  duration  'T ms  which  is  greater  than  TT  and 
is  given  by 
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where 
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This  extended  duration  is  then  used  in  Equation  (5.10)  to  estimate  N  and 
in  Equation  (5.13)  for  the  RMS  calculation. 

In  order  to  accommodate  surface  wave  effects  at  large  distances  (R>100  km) 
the  code  has  an  option  to  incorporate  R-^  geometrical  attenuation  rather 
than  the  body  wave  R-1  attenuation  beyond  100  km. 

The  dominance  of  higher  mode  surface  waves  beyond  about  two  crustal 
thicknesses  (JflOO  km)  also  cause  the  ground  motion  duration  to  be  longer 


than  the  source  duration  (  T  *  £c  ).  A  correction  factor,  based  upon 
empirical  and  synthetic  records,  has  been  developed  (Herrmann,  1985)  which 
shows  the  duration  proportional  to  distance 


'Y  ’  f'1  +  0  o  S'  R  (k^) 

(5.19) 


Both  of  these  distance  correction  should  be  used  for  distances  greater 
than  100  km. 

In  order  to  curtail  the  Brune  spectrum  at  the  high  frequencies,  a  high  cut 
filter  that  accounts  for  the  observation  that  acceleration  spectra  often 
show  a  sharp  decrease  with  increasing  frequency,  above  some  cut-off 
frequency  (fmax)  is  added  to  Equation  (5.1).  It  represents  a  near  site 
(Hanks,  1  982)  or  source  effect  (Papagerogiou  and  Aki ,  1  983)  which  curtails 
the  high  frequency  energy.  Its  effect  is  generally  modeled  as  a  low  pass 
fourth  order  Butterworth  filter  (Boore,  1983)  with  a  modulus  given  by 
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An  alternate  model  for  this  high  frequency  attenuation  is  of  exponential 
form. 


Exp  (-^fk) 


(Anderson  and  Hough,  1984), 


where  kappa  (k)  is  a  region  specific  parameter.  Incorporation  of  this 
filter  into  the  RVT  source  model  has  not  demonstrated  as  optimum  a  fit  to 


observed  peak  acceleration  data  as  the  fmax  type  filter  (Luco,  1985; 

Boore,  1986). 

In  many  ground  response  analyses,  there  is  an  interest  in  ground  motion  at 
embedment  depth  rather  than  at  the  free  surface.  In  order  to  accommodate 
this,  the  transfer  function  for  a  single  layerjassuming  normal  incidence, 
of  thickness  H,  has  been  incorporated  into  the  code.  All  calculations; 
peak  values,  time  histories,  and  response  spectra,  are  then  evaluated  at  a 
depth  h  (may  be  zero)  within  the  layer.  The  layer  is  defined  by  a  shear 
velocity,  density,  and  damping  and  is  assumed  to  overlie  an  elastic  half¬ 
space.  While  this  is  not  exact  since  the  half-space  Q  is  not  infinite  (a 
WUS  Q  of  300  is  probably  a  minimum  value)  the  error  should  not  exceed  a 
few  percent. 

The  transfer  function  for  the  layer  is  given  by 


C^k‘W 


\A>  B, 


where 

k.  .  f  ( i  ♦  -  f) 

=  layer  shear  wave  velocity 
^  =  dampi ng  of  1 ayer 

H 


(5.21) 


layer  thickness 


and  subscripts  1  and  2  refer  to  the  layer  and  half-space  respectively. 

For  non-zero  h  ,  only  the  modulus  of  HL  is  used  since  the  RVT  technique 
employs  the  power  spectrum  to  calculate  the  RMS  values  (Equation  (5.6)). 
For  time  domain  calculations,  this  may  introduce  some  signal  distortion 
but  generally  site  thicknesses  are  less  than  the  wavelengths  of 
predominant  energy  so  the  effects  should  be  minimal.  Several  tests  have 
shown  this  to  be  the  case. 

An  additional  caution  in  using  the  site  response  in  the  RVT  calculations 
arises  regarding  the  duration  calculations  (Herrmann,  1985).  If  a  very 
soft  site  is  specified,  such  that  considerable  energy  can  be  trapped  due 
to  reflected  waves,  the  effective  duration  is  increased  much  like  the 
higher  mode  surface  waves  at  large  distances.  This  effect  however,  is  not 
compensated  for  in  the  RVT  theory  so  the  peak  values  may  be  biased.  The 
RVT  approach  then  should  not  be  employed  to  model  expected  peak  or 
response  spectral  values  for  application's  where  site  response  dominates 
the  motion. 

Another  correction  factor,  designed  to  explicitly  account  for  wave 
ampl ification  due  to  decreasing  seismic  velocities  near  the  earth's 
surface,  has  been  incorporated.  Boore  (1985)  has  calculated  typical  WUS 
frequency  dependent  factors  based  upon  an  average  velocity  change  over  one 
quarter  wavelength  (Joyner  and  Fumal,  1984).  The  equation,  based  upon 
energy  conservation  is  given  by 


yv 


where  the  subscripts  o  and  R  refer  to  average  crustal  properties  and  near 
the  receiver  respectively.  The  WUS  factors  are  given  in  Table  2  for 
several  frequencies  (Boore,  1985).  Values  for  intermediate  frequencies 
are  linearly  interpolated  while  values  for  frequencies  which  are  outside 
the  specified  range  are  merely  extended. 


7.0  USER  MANUAL 


The  code  has  been  given  the  name  RASCAL  which  stands  for  Response  Spectra 
and  Acceleration  Scaling.  We  have  referred  to  it  by  many  less  eloquent 
names  during  its  development  but  RASCAL  seems  to  capture  the  essential 
elements.  No  plotting  routines  are  incorporated  since  these  are  highly 
user  and  machine  dependent.  Pertinent  files  are  created  and  the  user  is 
free  to  plot  sizes  and  the  scales  of  choice. 

Included  with  the  code  is  a  library  of  basis  time  histories  which  are  used 
to  generate  the  phase  spectra.  The  bases  are  shown  in  Figure  set  12  and 
each  plot  title  describes  the  time  history.  The  distance  and  magnitude 
selection  criterion  are  shown  in  Table  4.  The  basis  library  time  histo¬ 
ries  are  coded  in  local  magnitude  (M^)  and  in  hypocentral  range  (R).  The 
code  is  presently  compiled  in  FORTRAN  F77  on  a  PRIME  750  and  an  AT  compat¬ 
ible  microcomputer.  It  was  designed  to  be  as  machine  independent  as 
possible  to  allow  for  greatest  flexibility. 

7.1  Prediction  Methodology 

To  predict  peak  values  and  RVT  response  spectra,  relevant  parameters  must 
be  specified  along  with  the  proper  KEYS.  This  is  explained  in  detail  in 
the  User's  Guide  (7.4). 

7.2  Synthesis  Methodology 

Once  the  source  and  wave  propagation  parameters  are  specifed,  selection  of 
the  proper  KEY  (see  User's  Guide)  will  produce  a  time  history.  The  basis 
time  history,  from  which  the  phase  is  extracted,  may  be  user  supplied 
(external)  or  automatically  selected  from  the  library  (internal).  Table  4 
lists  the  library  contents  along  with  the  magnitude  and  distance  selection 
criterion.  If  a  time  history  is  supplied,  the  total  number  of  points  must 
not  exceed  2048.  Also,  since  both  an  RVT  response  spectrum  and  one  based 
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upon  a  time  domain  integration  are  produced,  the  total  length  of  the  time 
history  must  exceed  10  seconds.  This  arises  because  the  longest  period 
oscillator  in  the  time  domain  response  spectrum  calculation  is  10 
seconds.  Zeros  may  be  appended  prior  to  inputting  the  external  basis  to 
fill  out  the  array  to  10  seconds.  At  the  high  frequency  end,  the  highest 
oscillator  frequency  is  set  to  0.7  fn  (fn  =  Nyquest  frequency)  or  34  Hz, 
whichever  is  smaller. 

The  actual  oscillator  frequencies  are  included  in  a  data  statement.  They 
consist  of  the  NRC  specified  values  in  additional  to  selected  intermediate 
values  and  represent  a  standard  WCC  (Woodward-Clyde  Consultants)  format  of 
143  values  (see  User's  Guide). 

The  peak  time  domain  values  of  the  acceleration  or  velocity  are  normalized 
to  either  the  respective  RVT  predicted  peaks  or  to  input  values  (see 
User's  Guide).  When  integrating  to  velocity  or  displacement,  filtering 
should  be  performed.  The  filter  supplied  is  a  causal  Butterworth  which 
can  be  applied  twice  to  provide  band-pass  capability.  Guidelines  for 
corner  periods  for  the  high-pass  flter  vary  but  we  generally  assume  a 
sufficient  signal  to  noise  ratio  out  to  the  source  corner  period  (fc_1). 

We  therefore  put  the  filter  corner  at  a  slightly  longer  period  and  use  a 
fall-off  of  30  db/octave  (5  pole). 

7.3  Response  Spectral  Scaling  Methodology 

Either  pseudo  acceleration  or  psuedo  velocity  may  be  entered  as  the  target 
response  spectrum.  The  frequency  range  should  be  from  0.1  Hz  to  34  Hz 
with  enough  points  to  define  the  curve  (max  =  143).  Values  higher  than  34 
Hz  may  be  input,  but  only  values  between  0.1  Hz  and  34  Hz  are  employed. 

The  low  frequency  limit  (0.1  Hz)  is  fixed. 

The  first  few  (user  specified)  iterations  are  performed  using  the  RVT 
response  spectrum  to  scale  the  Brune  spectrum.  The  final  few  (user 
specified)  iterations  use  a  standard  time  domain  integration  algorithm  to 
calculate  the  response  spectrum  from  the  synthesized  acceleration  time 


history.  This  response  spectrum  is  then  used  to  scale  the  Brune 
modulus.  The  final  acceleration  time  history  may  be  either  normalized  to 
an  input  value  or  left  unsealed. 


PROGRAM  RASCAL 


PROGRAM  RASCAL  (RESPONSE  SPECTRUM  AND  ACCELEROGRAM  SCALINO) 
COMPUTES  BRUNE  FOURIER  AMPLITUDE  BASED  ON  THE  METHOD  SUGGESTED 
BY  J.  N.  BRUNE  (B38A,  VOL.  79.  SEPTEMBER  1970)  AND  COMPUTES  PEAK 
ACCELERATION.  PEAK  VELOCITY.  RESPONSE  SPECTRAL  ACCELERATION  AND 
RESPONSE  SPECTRAL  VELOCITY  BASED  ON  THE  METHOD  SUOOESTED  BY  D.  M. 
BOORE  (B88A.  VOL.  73.  DECEMBER  1983)  BY  USING  RANDOM  VIBRATION 
THEORY  (RVT)  TECHNIQUES. 

THE  PROORAM  AL80  GENERATES  SYNTHETIC  TIME  HI8T0RY  (ACCELERATION. 
VELOCITY  OR  DISPLACEMENT)  BY  COMPUTING  THE  FFT  AND  EXTRACTING  THE 
PHASE  OF  AN  INPUT  ACCELEROGRAM  (INTERNAL  OR  EXTREMAL  BASE)  AND 
COMBINES  THE  COMPUTED  BRUNE  FOURIER  AMPLITUDE  TO  GENERATE  THE 
OUTPUT  TIME  HISTORIES. 

FOR  RESPONSE  8PECTRUM  AND  ACCELEROORAM  SCALINO,  THE  PROORAM  SCALES 
THE  COMPUTED  RVT  RESPONSE  SPECTRUM  WITH  THE  INPUT  DESION  (TAROET) 
RESPONSE  SPECTRUM.  THIS  PROCESS  REPEATS  FOR  A  SPECIFIIED  NUMBER  OF 
TIMES  AND  THE  SPECTRAL  MATCHINO  SHOULD  CQNVERQE  AFTER  2  TO  3 
ITERATIONS.  DURING  THIS  ITERATIVE  SPECTRAL  SCALINO  OF  THE  RVT 
RESPONSE  SPECTRUM.  THE  COMPUTED  BRUNE  FOURIER  AMPLITUDE  IS  ALSO 
SCALED.  THE  PROORAM  THEN  COMBINES  THIS  SCALED  BRUNE  FOURIER 
AMPLITUDE  WITH  THE  PHASE  OF  THE  INPUT  ACCELEROORAM  AND  COMPUTES  AN 
OUTPUT  ACCELEROORAM.  THIS  OUTPUT  ACCELEROORAM  IS  THEN  USED  AS  AN 
INPUT  ACCELEROORAM  TO  COMPUTE  A  S I NOLE— DEOR EE -OF -FREEDOM  (SDF ) 
RESPONSE  SPECTRUM.  AGAIN,  THE  PROORAM  SCALES  THE  COMPUTED  SDF 
RESPONSE  SPECTRUM  WITH  THE  INPUT  (TAROET)  RESPONSE  SPECTRUM. 

SCALES  THE  BRUNE  FOURIER  AMPLITUDE  SIMULTANEOUSLY,  AND  COMPUTES  AN 
OUTPUT  ACCELEROORAM.  THIS  PROCESS  ALSO  REPEATS  FOR  ANOTHER 
SPECIFIED  NUMBER  OF  TIMES.  FINALLY,  THE  PROORAM  NORMALIZES  OR 
FILTERS  THE  OUTPUT  ACCELROORAM  AND  COMPUTES  SDF  RESPONSE  SPECTRUM 
BASED  ON  THIS  CONDITIONED  OUTPUT  ACCELEROORAM. 

THIS  PROORAM  HAS  THE  CAPABILITY  OF  COMPUTING  THE  HALF-SPACE  (ROCK 
OUTCROP)  OR  SITE  ( SINGLE-LAYERED  SOIL)  RESPONSE  AND  OF  GENERATING 
OUTPUT  VELOCITY  OR  DISPLACEMENT  TIME  HISTORY  AS  WELL. 


USER'S  GUIDE 
BASIC  PARAMETERS  : 

1  TITLE  :  FORMAT (A80) 

TITLE  IS  THE  TITLE  OF  THE  RUN. 

2  OFILE  :  F0MAT(A80> 

OFILE  IS  THE  OUTPUT  FILE  NAME  OF  THE  RUN 

OUTPUT  FILES  OF  SPECTRAL  VALUES,  ACCELERATION,  VELOCITY, 
AND  DISPLACEMENT  TIME  HISTORIES  ARE  AUTOMATICALLY  CREATED 
SEPARATELY  UNDER  THIS  OFILE  NAME  WITH  SUFFIXES  OF  BP,  A, 

V  AND  D  RESPECTIVELY.  AND  OF  NUMBER  OF  ITERATIONS  AT  THE 
END  OF  EACH  FILE. 

3  SDR OP.  DENS.  D.  H,  SV  :  FORMAT (FREE ) 
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8 DR OP  18  THE  BTRES8  DROP  (BARS) 

DENS  18  THE  MEDIUM  DENSITY  (OM/CC ) 

D  18  THE  EP I CENTRAL  DISTANCE  (KM). 

H  IS  THE  SOURCE  DEPTH  (KM). 

BV  IS  THE  SHEAR  WAVE  VELOCITY  OP  THE  HALF-SPACE  (KM/SEC). 

4.  Q.  FQ.  ALPHA.  MW  :  FORMAT  (FREE) 

Q  IS  THE  FREQUENCY  DEPENDENT  QUALITY  FACTOR  OF  THE  HALF-SPACE 
FOR  Q  FILTER.  WHERE  G(F)  -  Q  *  (  F  /  FQ  )  **  ALPHA. 

FQ  IS  THE  CONTROL  FREQUENCY  FOR  THE  Q  QUALITY  FACTOR  (HERTZ). 
ALPHA  IS  THE  EXPONENTIAL  CONSTANT  FOR  THE  Q  OPERATOR. 

MW  IS  THE  MOMENT  MAGNITUDE. 

9.  DAMP.  FMAX.  N.  CAP  :  FORMAT  (FREE) 

DAMP  IS  THE  SPECTRAL  DATING  (PERCENT). 

(MU8T  BE  BETWEEN  1  TO  10  X). 

FMAX  IS  THE  BUTTERWORTH  CORNER  FREQUENCY  (HERTZ). 

(IF  FMAX  -  0.  NO  FILTERING  18  PERFORMED). 

N  IS  THE  FMAX  ORDER  NUMBER. 

(IF  N  -  0.  NO  FILTERING  IS  PERFORMED). 

CAP  IS  THE  CONSTANT.  KAPPA  FOR  NEAR-SITE  EXPONENTIAL 
FILTERING. 

(IF  CAP  -00,  NO  FILTERING  IS  PERFORMED). 

FUNCTIONAL  KEYS  : 

6.  KEY1 ,  KEY2.  KEY3,  KEY4.  KEY9.  KEY6  :  FORMAT (FREE) 

KEY1  IS  THE  KEY  FOR  COMPUTING  SPECIFIED  RESPONSE. 

IF  KEY1  -  O  :  COMPUTES  HALF-SPACE  RESPONSE. 

IF  KEY1  -  1  :  COMPUTES  SITE  RESPONSE. 

KEY2  IS  THE  KEY  FOR  APPLYINO  NEAR-SITE  AMPLIFICATION  FACTOR 
AS  A  FUNCTION  OF  FREQUENCY. 

IF  KEY2  -  O  :  NO  AMPLIFICATION  FACTORS  ARE  APPLIED. 

IF  KEY2  -  1  :  AMPLIFICATION  FACTORS  ARE  APPLIED 
KEY3  IS  THE  KEY  FOR  CHANGING  GEOMETRICAL  ATTENUATION  FOR 
R  >  lOO  KM 

IF  KEY3  -  O  BODY  WAVE  (  1  /  R  >  IS  RETAINED  FOR 
FOR  R  >  100  KM. 

IF  KEY3  -  1  :  ATTENUATION  FOR  R  >  100  KM  IS  APPLIED. 

(  FOR  R  <  100  KM.  .  ATTENUATION  -  1  /  R  > 

(  FOR  R  >  100  KM.  ,  ATTENUATION  -  1  /  100  * 
SORT (  100  /  R  )  ) . 

KEY4  IS  THE  KEY  FOR  INCREASING  EFFECTIVE  SOURCE  DURATION  DUE 
TO  SURFACE  WAVE  CONTRIBUTION. 

IF  KEY 4  -  0  DURATION  -  1  /  FC  . 

IF  KEY4  -  1  DURATION  -  1  /  FC  ♦  0.  09  *  R  . 

(FROM  R  B  HERRMANN,  BSSA,  VOL.  79.  OCTOBER  1989) 
KEY9  IS  THE  KEY  FOR  OENERATINO  ARTIFICIAL  ACCELERATION  TIME 
HISTORY 

IF  KEYS  -  0  NO  ARTIFICIAL  ACCELERATION  TIME  HISTORY  IS 
GENERATED  JUST  BASIC  RESPONSE  IS  COMPUTED 
IN  THIS  CASE,  NO  NUMBER  OF  ITERATIONS  AND 
DESION  (TAROET)  RESPONSE  SPECTRUM  ARE 
REQUIRED.  SET  KEY6  -  0  ! ! ! 

IF  KEYS  >  l  :  ARTIFICIAL  ACCELERATION  TIME  HISTORY  IS 
GENERATED  USING  PHASES  OF  THE  BUILT-IN 
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ACCELERATION  TIME  HI9T0RY. 

IF  KEYS  <  1  :  ARTIFICIAL  ACCELERATION  TIME  HI8T0RY  IS 
GENERATED  USING  PHASE8  OF  THE  INPUT 
ACCELERATION  TIME  HISTORY. 

KEY 6  IS  THE  KEY  FOR  WRITING  OUTPUT  TIME  HISTORIES.  OFILE 

IF  KEY*  -  0  :  ONLY  OUTPUT  ACCELERATION  TIME  HISTORY  OF  THE 
LAST  ITERATION  IS  WRITTEN  IN  OFILE.  A.  OR  NO 
OUTPUT  TIME  HISTORIES  ARE  WRITTEN  AT  ALL  IF 
KEYS  -  0. 

IF  KEY*  -  1  :  ONLY  OUTPUT  ACCELERATION  AND  VELOCITY  TIME 
HISTORIES  ARE  WRITTEN  IN  OFILE.  A  AND 
OFILE.  V.  RESPECTIVELY. 

IF  KEY*  -  2  :  ONLY  OUTPUT  ACCELERATION  AND  DISPLACEMENT 
TIME  HISTORIES  ARE  WRITTEN  IN  OFILE. A  AND 
OFILE.  D,  RESPECTIVELY. 

IF  KEY 6  -  3  :  ALL  OUTPUT  ACCELERATION.  VELOCITY  AND 

DISPLACEMENT  TIME  HISTORIES  ARE  WRITTEN  IN 
OFILE.  A.  OFILE.  V  AND  OFILE.  D.  RESPECTIVELY 

PARAMETERS  FOR  THE  FUNCTION  KEYS  : 

SOIL  RESPONSE  :  SKIP  THIS  GROUP  7  IF  KEY1  -  0  ! ! ! 

7.  1  DENI.  SV1.  DAMP  1 .  THICK.  RD  :  FORMAT (FREE) 

DENI  IS  THE  SOIL  DENSITY  (CM/CC). 

SV1  IS  THE  SHEAR  WAVE  VELOCITY  OF  THE  SOIL  (M/SEC). 

DAMP1  IS  THE  DAMPINO  RATIO  OF  THE  SOIL  (PERCENT) 

THICK  IS  THE  THICKNESS  OF  THE  SOIL  LAYER  (M) 

RD  IS  THE  RECEIVER  DEPTH  (DEPTH  OF  INTEREST)  (M) 

(RD  MUST  BE  LE.  THICK). 

72  DEN2.  SV2  :  FORMAT (FREE) 

DEN2  IS  THE  BEDROCK  DENSITY  (GM/CC) 

SV2  IS  THE  SHEAR  WAVE  VELOCITY  OF  THE  BEDROCK  (M/SEC) 

***  NOTE  :  DENSITY  UNITS  IN  DENI  AND  DEN2  CAN  BE  IN  OTHER 
UNITS  IF  CONSISTENT  ! ! ! 

LENGTH  UNITS  IN  THICK.  RD.  SV1  AND  SV2  CAN  BE  IN 
OTHER  UNITS  IF  CONSISTENT  ! ! ! 

NEAR-SITE  AMPLIFICATION  FACTORS  :  SKIP  THIS  GROUP  8  IF  KEY2  =*  O  !  !  ! 

0  NAMP.  AFREQ(I),  AMFAC ( I )  :  FORMAT (FREE) 

NAMP  IS  THE  NUMBER  OF  AMPLIFICATION  FACTORS 

AFREO(I)  IS  THE  FREQUENCY  AT  WHICH  AMPLIFICATION  IS  APPLIED 
(HERTZ) 

AMFAC ( I )  IS  THE  AMPLIFICATION  FACTOR  AT  AFREQ(I) 

***  REPEAT  (AFREQ(I). AMFAC ( I ) )  PAIR  NAMP  TIMES.  (MAXIMUM  100)  ! 

***  NOTE  .  FOR  FREQUENCY  <  AFREQ(l),  AMFAC ( 1 )  IS  USED 

FOR  FREQUENCY  >  AFREQ ( NAMP ) .  AMFAC ( NAMP )  IS  USED 

INPUT  ACCELERATION  TIME  HISTORY  SKIP  THIS  OROUP  9  IF  KEYS  GE  0 

9  1  AFILE  FORMAT ( A80 ) 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


AFILE  IS  THE  INPUT  ACCELERATION  TIME  HI8T0RY  FILE  NAME. 

9  2  NPT,  NHEAD.  DT,  FMT  :  F0RMAT(2I3,  FIO.  4.  A40> 

NPT  18  THE  NUMBER  OF  PGINT8  IN  THE  INPUT  ACCELERATION  TIME 
HISTORY.  (MAXIMUM  4096). 

NHEAD  18  THE  NUMBER  OF  HEADER  CARDS  IN  THE  INPUT  ACCELERATION 
TIME  HISTORY. 

DT  IS  THE  TIME  INCREMENT  OF  THE  INPUT  ACCELERATION  TIME 
HISTORY  (SECONDS). 

FMT  IS  THE  READ  FORMAT  OF  THE  INPUT  ACCELERATION  TIME  HISTORY. 
DEFAULT  IS  (8F9.  6). 

***  NOTE  :  IF  NPT  <  2048.  TRAILING  ZEROS  ARE  AUTOMATICALLY 
ADDED  TO  DOUBLE  THE  NPT  TO  THE  CLOSEST  POWER  OF 
2  POINTS. 

IF  NPT  >  2048.  ONLY  THE  FIRST  2048  POINTS  ARE  USED. 
AND  2048  POINTS  OF  TRAILING  ZER08  ARE  AUTOMATICALLY 
ADDED  TO  THE  TOTAL  OF  4096  POINTS. 

-  OUTPUT  TIME  HISTORIES  :  SKIP  THIS  CROUP  10  IF  KEY3  -  O  ! ! ! 

10.  1  FACTA.  FACTV.  FACTD  :  FORMAT (FREE) 

FACTA.  FACTV  fc  FACTD  ARE  THE  NORMALIZING  FACTORS  FOR  THE 
OUTPUT  ACCELERATION.  VELOCITY  fc  DISPLACEMENT  TIME 
HISTORIES.  RESPECTIVELY. 

IF  FACTA.  V.  D  -0.0:  NO  NORMALIZATION  IS  APPLIED. 

IF  FACTA.  V.  D  <0.0:  THE  OUTPUT  TIME  HISTORIES  ARE  NORMALIZED 

TO  THE  COMPUTED  PEAK  VALUES  BASED  ON  THE 
RANDOM  VIBRATION  THEORY.  RESPECTIVELY. 

IF  FACTA. V, D  >00:  THE  OUTPUT  TIME  HISTORIES  ARE  NORMALIZED 

TO  FACTA.  V.  D.  RESPECTIVELY. 

10  2  NIT1,  NIT2  :  FORMAT ( FREE ) 

NIT1  IS  THE  NUMBER  OF  ITERATIONS  USED  FOR  RVT  SPECTRAL  MATCHING 
NIT2  IS  THE  NUMBER  OF  ITERATIONS  USED  FOR  SDF  SPECTRAL  MATCHING 

***  NOTE  :  2  TO  3  ITERATIONS  FOR  EACH  CASE  IS  RECOMMENDED 

-  DESIGN  (TARGET)  RESPONSE  SPECTRUM  :  SKIP  THIS  QROUP  10.3 

IF  NIT1  AND  NT2  -  O  ! ( ! 

10  3  NPSV.  TFREQ(I).  TSV(I)  :  FORMAT (FREE) 

NPSV  IS  THE  NUMBER  OF  DESIGN  (TAROET)  SPECTRAL  VALUES 
(MAXIMUM  1O0) 

IF  NPSV  >  0  SPECTRAL  VALUES  (TSV(I))  MUST  BE  IN 
ACCELERATION  (G'S) 

IF  NPSV  <  O  SPECTRAL  VALUES  (TSV(I)>  MUST  BE  IN 
VELOCITY  (CM/SEC) 

TFREQ(I)  IS  THE  DESION  (TAROET)  FREQUENCY  (HERTZ). 

TSV(I)  IS  THE  DESIGN  (TARGET)  SPECTRAL  VALUES  AT  TFREQ(I) 

***  REPEAT  (TFREQ( I >, TSV < I ) )  PAIR  ABS(NPSV)  TIMES  M1 

***  NOTE  SET  THE  DESION  FREQUENCY  RANOE  BEYOND  O  1  TO  34  O 
HERTZ,  IF  POSSIBLE  I.E  SET  TFREQ  (  1  >  LIL  0.1 
AND  TFREQ ( NPSV >  QE  34  O  HERTZ.  THIS  RANOE  l* 


14 


UP  FOR  INTERPOLATING  THE  DESIGN  SPECTRAL  VALUES  AT 
THE  OSCILLATOR  FREQUENCY  REOIME  !'! 

-  FILTERING  PARAMETERS  FOR  COMPUTING  OUTPUT  TIME  HI8TQRIEB 

10.4  FC1.  NF1,  FC2,  NF2  :  FORMAT  (FREE) 

FC1  IS  THE  CORNER  FREQUENCY  OF  THE  FIRST  BUTTERWORTH  FILTER, 
(-VE  FOR  REMOVAL).  (HERTZ). 

NF1  IS  THE  ORDER  NUMBER  (NO.  OF  POLES)  OF  THE  FIRST 
BUTTERWORTH  FILTER.  (-VE  FOR  HIOH  PASS). 

FC2  IS  THE  CORNER  FREQUENCY  OF  THE  SECOND  BUTTERWORTH  FILTER, 
(-VE  FOR  REMOVAL).  (HERTZ). 

NF1  IS  THE  ORDER  NUMBER  (NO.  OF  POLES)  OF  THE  SECOND 
BUTTERWORTH  FILTER.  (-VE  FOR  HIOH  PASS). 
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-  PROGRAM  rascal 

PROGRAM  RASCAL  (RESPONSE  SPECTRUM  AND  ACCELEROGRAM  SCALING) 
COMPUTES  BRUNE  FOURIER  AMPLITUDE  BASED  ON  THE  METHOD  SUGGESTED 
BY  J.  N.  BRUNE  (B8SA,  VOL.  73,  SEPTEMBER  1970)  AND  COMPUTES  PEAK 
ACCELERATION,  PEAK  VELOCITY,  RESPONSE  SPECTRAL  ACCELERATION  AND 
RESPONSE  SPECTRAL  VELOCITY  BASED  ON  THE  METHOD  SUGGESTED  BY  D.  M 
BOORE  ( BSSA,  VOL.  73.  DECEMBER  1983)  BY  U8IN0  RANDOM  VIBRATION 
THEORY  (RVT)  TECHNIQUES 

THE  PROGRAM  ALSO  GENERATES  SYNTHETIC  TIME  HI8T0RY  (ACCELERATION, 
VELOCITY  OR  DISPLACEMENT)  BY  COMPUTING  THE  FFT  AND  EXTRACTING  THE 
PHASE  OF  AN  INPUT  ACCELEROGRAM  (INTERNAL  OR  EXTRENAL  BASE)  AND 
COMBINES  THE  COMPUTED  BRUNE  FOURIER  AMPLITUDE  TO  GENERATE  THE 
OUTPUT  TIME  HISTORIES. 

THIS  PROORAM  HAS  THE  CAPABILITY  OF  COMPUTING  THE  HALF-SPACE  (ROCK 
OUTCROP)  OR  SITE  (SINGLE-LAYERED  SOIL)  RESPONSE  AND  OF  GENERATING 
OUTPUT  VELOCITY  OR  DISPLACEMENT  TIME  HISTORY  AS  WELL. 

-  CODED  BY  KIN  W  LEE,  NOVEMBER,  1983. 

WOODWARD-CLYDE  CONSULTANTS.  WALNUT  CREEK  OFFICE 


CHARACTER *80 
St 
St 


TITLE.  OFILE,  AFILE,  WORD,  FLA0(143>*1,  FMT*40, 
FILES.  FILE9,  FILE11.  FILE12,  FILE13,  SUFB*4, 
SUF9*3.  SUFI 1*2.  SUFI 2*2,  SUFI 3*2,  C IT  1*2,  CIT*2 


REAL  MW 


DIMENSION  FREQ( 143 ) .  FAS(143>.  RSA(143).  RSV(143),  PAA(143), 

S,  PRV  ( 1 43  > .  TSVK145).  SPRAT(143),  TFREQ(IOO),  TSV(IOO) 


COMPLEX  RE.  HS 


COMMON  SFREQOOOl),  SFASOOOl).  FUNXOOOOl).  FUNX20001), 

&  FUNX4  <  3001 ) .  SWI <  3001 ) ,  SPRATS( 3001 ) , 

fc  CP ( 2030 ) .  SP ( 2050 ) ,  0TH(4100> 


COMMON 

/PI DATA/ 

PI,  P 12 

COMMON 

r 

/FADATA/ 

KEY1.  KEY2, 
CAP,  PICA P 

C.  FC.  PIR,  Q.  FQ. 

ALPHA, 

COMMON 

/AMDATA/ 

AFREQ( 100). 

AMFAC(IOO).  NAMP 

COMMON 

/SMDATA/ 

SV1.  DAMP  1 . 

THICK.  ARATIO,  RD, 

RK,  CK 

COMMON 

t 

/TMDATA/ 

NP I  A.  NPT .  NPT2.  NP2.  NPFR ,  DT, 

A  <  4 1 00 ) ,  PFREQ ( 2030 ) .  PFAS(20S0> 

COMMON 

/SPDATA/ 

DAMP.  NFREQ. 

WI(143).  PER (143) 

FMAX, 


N2, 


C 

C  -  OSCILLATOR  (WCC)  FREQUENCY  REOIME 
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DATA 

FREQ/O.  001. 

0.  0999, 

Sc 

0.  100. 

0.  111. 

0.  123. 

0.  143. 

0.  167, 

0.  182. 

0.  200, 

Sc 

0.  222. 

0.  230, 

0.  263, 

0.  278. 

0.  2R4. 

0.  300, 

0.  312, 

Sc 

0.  333. 

0.  337, 

0.  383. 

0.  400, 

0.  417, 

0.  433, 

0.  300. 

Sc 

0.  356. 

0.  600, 

0.  623. 

0  667, 

0.  700, 

0.  714. 

0  769, 

Sc 

0.  800. 

0  833, 

0.  900, 

0.  909, 

1.  OOO. 

1.  100, 

1.  Ill, 

Sc 

1.  176. 

1.  200, 

1.  230, 

1.  300, 

1.  333, 

1 .  400. 

1.  429, 

Sc 

1.  471. 

1.  300, 

1.  313. 

1  362, 

1.  600, 

1.  613, 

1  667, 

Sc 

1.  700. 

1.  724. 

1.  786. 

1  800, 

1.  832, 

1.  900, 

1  923. 

Sc 

2.  000. 

2.  083, 

2.  100, 

2.  174, 

2.  200, 

2.  273, 

2.  300, 

s. 

2.  381. 

2.  400, 

2.  300, 

2.  600, 

2.  632, 

2.  700. 

2  778, 

Sc 

2.  800. 

2  900, 

2.  941. 

3  OOO, 

3.  123. 

3.  133. 

3.  300, 

Sc 

3.  333. 

3.  448, 

3.  371, 

3.  600, 

3.  BOO. 

3.  830, 

4  OOO, 

Sc 

4.  167. 

4  200, 

4  400, 

4  330, 

4.  600. 

4.  800, 

3.  OOO, 

Sc 

3.  230. 

3.  263, 

3  300, 

3.  336. 

3.  730. 

3.  882. 

6.  OOO, 

Sc 

6.  230. 

6.  300, 

6.  667, 

6.  730, 

7.  OOO. 

7.  143, 

7.  230, 

Sc 

7  300, 

7  692, 

7  730, 

8  OOO, 

8.  333. 

8.  500, 

9.  OOO. 

Sc 

9  091. 

9  300, 

10  OOO, 

10.300,  11.000,  11 

111.  11.  300, 

Sc 

11.  763. 

12.  OOO 

i.  12.  300,  13. 

OOO,  13. 

333,  13.  300,  14.  OOO, 

Sc 

14.  286. 

14.  300 

>,  13  OOO,  13. 

383,  16 

OOO,  16  667,  17  OOO, 

Sc 

18.  000. 

18  868 

1.  20.  OOO,  22. 

OOO,  25 

OOO,  28  OOO,  29  412, 

Sc 

31  OOO. 

33.  333 

,  34.  000/ 

READ<5, 1100)  TITLE 

READO.  1100)  OFILE 

READO.  *>  SDROP,  DENS.  D.  H.  SV 

READO.*)  Q,  FQ,  ALPHA.  MW 

READO.  *)  DAMP,  FMAX,  N.  CAP 

READO,*)  KEY  1 ,  KEY2,  KEY3,  KEY4.  KEYS,  KEY6 

PRINT  1200,  TITLE 
PRINT  1300,  OFILE 

PRINT  1340,  SDROP,  DENS.  D,  H,  SV,  0,  FQ,  ALPHA,  MW,  DAMP, 

8.  FMAX,  N,  CAP,  KEY1 ,  KEY2,  KEY3,  KEY4,  KEY5,  KEY6 

FUNCTIONAL  KEY 1  OPERATION 

PRINT  1200.  TITLE 
PRINT  1400 

IF  (  KEY1  EQ  0  )  THEN 
PRINT  1410 
ELSE 

READO.*)  DENI.  SV1,  DAMP  1 ,  THICK,  RD 
READO.*)  DEN2,  SV2 

ARATIO  -  (  DENI  *  SV1  )  /  <  DEN2  *  SV2  ) 

PRINT  1420 

PRINT  1430.  DENI,  SV1,  DAMP  1 ,  THICK,  RD,  DEN2,  SV2 
DAMP  1  -  DAMP  1  /  100  0 
END  IF 

FUNCTIONAL  KEY2  OPERATION 
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oooo  noon  nnn  n  o  n  o 


■  '  i  w1  ail  ■] 


w 


C 


IF  (  KEY2  .  EQ.  0  )  THEN 
PRINT  1440 
ELSE 

READ  <5.  *>  NAMP,  <  AFREQCI),  AMFAC < I > .  1=1, NAMP  ) 

PRINT  1460 
DO  120.  1  =  1.  NAMP 

120  PRINT  1480.  I,  AFREQ < I ) .  AMFAC ( I ) 

END  IF 

-  INITIALIZES  CONSTANT  PARAMETERS 

NFREQ  -  145 

AMPF  -  1.0 

DAMP  »  DAMP  /  100.  0 

DAMP 2  »  DAMP  *  2  0 

N2  -  N  *  2 

PI  -  3.  141592653589793 

P 12  *  2  0  *  PI 

PICAP  =  PI  *  CAP 

SR  HP  I  -  SORT  (  0.5  *  PI  ) 

EULER  =  0  5772156649 

SDAMP  -  lOO  0  /  (  2.  0  *  Q  > 

FC  =  21.  07  *  SV  *  (  SDROP  **  (  1.0/30)  )  *  EXP(-1.  1 
SR  =  (  2.  34  *  SV  )  /  (  PI2  *  FC  ) 

R  -  SORT (  D  *  D  ♦  H  *  H  ) 

-  FUNCTIONAL  KEY3  OPERATION 

ATTEN  *  1  0  /  R 

IF  (  R  GT  lOO  O  AND  KEY3  .  EQ.  1  )  THEN 

ATTEN  =01/  SQRT(R) 

PRINT  1500,  ATTEN 
ELSE 

PRINT  1520.  ATTEN 
END  IF 

C  =  <  7  8  *  SDROP  *  SR  *  ATTEN  )  /  (  DENS  *  SV  ) 

-  FUNCTIONAL  KEY4  OPERATION 

HERR  =00 

IF  <  KEY4  EQ.  0  )  THEN 
PRINT  1540 
ELSE 

HERR  =  O  05  *  R 
PRINT  1560,  HERR 
END  IF 

DRVT  =  1  0  /  FC  +  HERR 
PRINT  1580,  DRVT 

-  FUNCTIONAL  KEYS  OPERATION 

NP2  =  12 
NPT2  =  4098 
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5  *  MW) 


r>  n  r>  n  nooo  o  o  r> 


C 


NPT  -  4096 
NHEAD  -  3 
DT  -  0.  01 
FMT  -  '<8F9.  6) ' 

C 

NPFR  -  2049 
PFR1  -  0.0244141 
DPFR  -  PFR1 
C 

FACTA  -  O.  O 
FAC  TV  -0.0 
FACTO  -  0.  0 
C 

TSVWX  -  O.  O 
DO  130  I-l.NFREQ 
130  T8VI  (  I  )  -0.0 
C 

IF  <  KEY 3  )  190.  140.  160 

-  KEY 3  -  O  :  NO  GENERATION  OF  SYNTHETIC  ACCELERATION  TIME  HISTORY 

140  PRINT  1600 
NIT1  -  O 
NIT  -  O 
CO  TO  340 

-  KEYS  >  O  :  GENERATION  OF  SYNTHETIC  ACCELERATION  TIME  HISTORY  BASED 

ON  A  BUILT-IN  ACCELERATION  TIME  HISTORY 

160  PRINT  1620.  KEY3 

-  CHOOSES  A  BUILT-IN  ACCELERATION  TIME  HISTDRY  BASED  ON  MW  AND  R 

RANGES 

IF  (  MW  .  LE.  4.  3  )  THEN 
IF  <  R  .  LE.  30.  0  )  THEN 
AFILE  -  'M4.NF.  DATA' 

NP I A  -  1320 
OO  TO  170 
ELSE 

AFILE  -  'M4.  FF.  DATA' 

NPIA  -  1320 
CO  TO  170 
END  IF 
END  IF 
C 

IF  (  MW  .  GT  4.  3  .  AND.  MW  LE  3.  3  )  THEN 
IF  (  R  .  LE.  30.  0  )  THEN 
AFILE  «  'M3.  NF.  DATA' 

NPIA  -  1320 
CO  TO  170 
ELSE 

AFILE  «  'M3.  FF  DATA' 

NPIA  -  2048 
GO  TO  170 
END  IF 
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c 


c 

c - 

c 

c 

190 

c 


END  IF 

IF  <  HU  .  OT.  3.  S  .  AND.  MW  .  LE.  6.  5  )  THEN 
IF  <  R  .  LE.  30.  0  )  THEN 
AFILE  -  'HA.  NF.  DATA' 

NPIA  -  2048 
00  TO  170 
ELSE 

AFILE  -  'M6.  FF.  DATA' 

NPIA  -  2048 
00  TO  170 
END  IF 
END  IF 

IF  <  R  .  LE.  30.  O  )  THEN 
AFILE  -  'H7.  NF.  DATA ' 

NPIA  -  2048 
ELSE 

AFILE  -  'M7.  FF.  DATA' 

NPIA  -  2048 
END  IF 

OPEN (7.  FILE-AFILE) 

DO  180  1-1, NHEAD 

READ <7, llOO)  WORD 

PRINT  1680,  WORD 

READI7.  FMT)  (  A<I>,  I-l.NPIA  ) 

CLOSE! 7) 

PRINT  1640,  AFILE.  MW.  R,  NPIA.  DT 
CO  TO  280 

KEYS  <  O  GENERATION  OF  SYNTHETIC  ACCELERATION  TIME  HISTORY  BASED 
ON  THE  INPUT  ACCELERATION  TIME  HISTORY 

PRINT  1660.  KEYS 


READO.  llOO)  AFILE 

READCS. 1700)  NPIA.  NHEAD,  DT,  FMT 

PRINT  1640.  AFILE.  MW.  R,  NPIA.  DT 

IF  (  NPIA  OT  2048  >  THEN 
NPIA  -  2048 
PRINT  1710 
END  IF 


IF  (  FMT  EQ 


)  FMT 


'  ( BF9.  6) 


C 

C - 

C 


OPEN ( 7 ,  FILE-AFILE) 

DO  220  I >1.  NHEAD 

READ (7. 1100)  WORD 

PRINT  1680.  WORD 

READ (7.  FMT)  (  A(I).  I-l.NPIA) 

CLOSE! 7) 

REDEFINE  NPT  IF  NPIA  <  2048 


NP2  -  HINT (  ALOO ( FLOAT < NP I A ) >  /  AL00(2.  0)  O  .  3  > 

NP2  -  NP2  ♦  1 
IF  (  NP2  .  CT.  12  )  THEN 
PRINT  1720,  NP2 
NP2  -  12 
END  IF 

NPT  -  2  **  NP2 
NPT2  -  NPT  ♦  2 

NPFR  -  NPT  /  2  ♦  1 
PFR1  -  1.  0  /  (  NPT  *  DT  ) 

DPFR  -  PFR1 

-  DETERMINES  THE  NYQUIST  FREQUENCY  OF  THE  INPUT  ACCELERATION  TIME 

HISTORY  AND  THE  UPPER  LIMIT  OF  THE  RESPONSE  SPECTRAL  FREQUENCY 
BY  TAKING  70  X  OF  THE  NYQUIST  FREQUENCY. 

FNYQ  -  O.  3  /  DT 
FNYO  -  O.  7  *  FNYQ 
IF  <  FNYQ  .  LT.  34.  O  )  THEN 
DO  240  K-l.NFREQ 

IF  (  FNYQ  .  LT.  FREQ(K)  )  00  TO  260 
240  CONTINUE 
260  NFREQ  -  K  -  1 
END  IF 

280  PRINT  1730.  NPFR.  PFR1.  DPFR 

-  READS  IN  NORMALIZING  FACTORS,  OR  DESIGN  (TARGET)  RESPONSE  SPECTRUM 

READO.  •>  FACTA.  FAC  TV.  FACTD 
PRINT  1740.  FACTA.  FAC TV.  FACTD 

READO.  •>  NIT!.  NIT2 

NIT  »  NIT1  ♦  NIT2 

PRINT  1700.  NIT.  NIT1.  NIT2 

IF  (  NIT  EQ.  0  )  GO  TO  340 

PRINT  1200.  TITLE 
PRINT  1760 

READO.  •>  NPSV.  (  TFREQ(I).  TSV(I),  1-1 ,  ABS(NPSV)  > 

IF  <  NPSV  OT  0  )  THEN 
PRINT  1770 
ELSE 

PRINT  1780 
END  IF 

DO  300  1-1.  ABS(NPSV) 

300  PRINT  1480.  I.  TFREQ(I).  TSV<I> 

-  INTERPOLATES  DESIGN  (TARGET)  RESPONSE  SPECTRUM.  TSV  AT  THE 

OSCILLATOR  FREQUENCY,  FREQ(I) 


non  non 


C 


DO  320  1-1,  NFREQ 

CALL  INTERP(  TFREQ.  T8V.  ABS(NPSV).  FREQ ( I  > »  BSV  ) 

TSVKI)  -  BSV 

IF  <  TSVHX  .  LE.  TSVKI)  )  TSVMX  -  TSVI(I) 

320  CONTINUE 

-  INITIALIZES  NUMBER  OF  ITERATIONS  TO  SUFFIXES,  CIT1  AND  CIT 

340  0PEN(21.  STATUS- 'SCRATCH' ) 

WRITE(21, '(212)')  NIT1 ,  NIT 
REWIND  21 

READ (21,  '(2A2)')  CIT1.  CIT 
IF  <  CITK1.1)  .EQ.  '  '  )  CIT1  ( 1:1)  -  'O' 

IF  <  CIT(1:  1)  .  EQ.  '  '  >  CIT(l.l)  -  'O' 

CLOSE (21 > 

-  FUNCTIONAL  KEY6  OPERATION 

I  -  INDEX (OFILE.  '  ') 

C 

K  -  1 

BUFB  -  RVT '  • 

FILES  -  OFILE( 1: I— 1 ) //SUF8//C IT1 
0PEN(8.  FILE-FILEB) 

C 

IF  (  KEYS  .  EQ.  O  )  THEN 
PRINT  1790,  FILE8 
00  TO  400 
END  IF 

C 

PRINT  1800.  KEY6 
PRINT  1820.  K.  FILES 
C 

K  *  K  +  1 
SUF9  -  'SP ' 

FILE9  -  0FILE(1: 1-1 ) //SUF9//CIT 
PRINT  1820.  K.  FILE9 
QPEN(9.  FILE-FILE9) 

C 

K  -  K  ♦  1 
SUFI  1  -  'A' 

FILE11  *  0FILE(1.  I-D//SUF11//CIT 
PRINT  1820.  K.  FILEU 
0PEN(11.  FILE-FILE1 1 ) 

C 

IF  (  KEY6  EQ.  0  )  GO  TO  390 
C 

CO  TO  <  360.  380.  360  )  KEY6 
C 

360  K  -  K  ♦  1 

SUFI 2  *  'V' 

FI  LEI  2  -  OFILEd  :  1-1  ) //SUF12//C  IT 
PRINT  1820.  K.  FILE 12 
OPEN( 12.  FILE-FILE12) 

IF  (  KEY6  .  NE  3  )  GO  TO  390 
C 
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<■>  o  r>  r>  n  onooooo  o  oooo  o  ooo 


360  K  -  K  ♦  1 

SUFI 3  -  D' 

FILE13  -  0FILE(1: 1-1 >//SUF13//CIT 
PRINT  1820.  K.  FI LEI 3 
0PEN( 13.  FILE-FILE13) 

390  READ (3.  *)  FC1.  NF1.  FC2.  NF2 
PRINT  1830.  FC1.  NF1,  FC2.  NF2 


400  CONTINUE 

PRINT  1200.  TITLE 

PRINT  1840.  SDAMP,  FC.  SR.  R 

PIR  -  PI  *  R  /  SV 
C  -  C  /  981.0 

-  INITIALIZES  RELATED  ARRAYS.  BASICALLY  FOR  PEAK  ACCELERATION  AND 

VELOCITY  COMPUTATIONS 

DO  420  1-1.  NFREQ 
FAS(I)  -0.0 
RSA(I)  -  0.  0  . 

RSV(I)  -00 
SPRAT  (  I  )  -1.0 
PER ( I  )  -0.0 
FLAG ( I )  -  '  ' 

420  CONTINUE 

DO  440  1-3.  NFREQ 

PER <  I  )  -1.0/  FREQ(I) 

WI(I)  -  PI2  *  FREO(I) 

440  CONTINUE 

-  COMPUTES  FOURIER  AMPLITUDE.  FAS  IN  THE  OSCILLATOR  FREQUENCY  REOIME 

CALL  BRUNEI  3.  NFREQ.  FREQ.  FAS  > 

-  COMPUTES  FOURIER  AMPLITUDE,  SFAS  IN  THE  BRUNE'S  OR  INTEGRATION 

FREQUENCY  REOIME 

NSFR  -  3001 
SFR1  -  O.  OOOl 
DSFR  -  O.  01 

DO  300  I  —  1 ,  NSFR 

SFREQ(I)  -  SFR 1  +  DSFR  *  FLOAT (  I  -  1  ) 

SWIII)  -  PI2  *  SFREQ(I) 

SPRATS!  I)  -1.0 
300  CONTINUE 

CALL  BRUNEI  1.  NSFR,  SFREQ,  SFAS  ) 


-  COMPUTES  RESPONSE  SPECTRAL  ACCELERATION,  RSA  AND  VELOCITY,  RSV 


00  non  n  non  nnnn  n  non  n  non  o  n n n  nn 


BY  RANDOM  VIBRATION  THEORY  <RVT>  FOR  NITI  NUMBER  OF  ITERATIONS 
PRINT  1860 

-  INITIALIZES  MAXIMUM  VALUES 

FASMX  -0.0 
RSAMX  -  0.  0 
RSVMX  -0.0 
SPRMX  -  0.  0 

DO  860  IT-1.  NITI  +  1 

-  UPDATES  8FAS(I)  WITH  NEW  SPRATS ( I > 

DO  620  I— 1. NSFR 

8FAS( I )  -  SFAS(I)  •  SPRATS ( I > 

620  CONTINUE 

DO  800  I-l.NFREO 

-  UPDATES  FAS ( I )  WITH  NEW  SPRAT ( I ) 

FAS(I)  -  FAS< I >  *  SPRAT ( I ) 

IF  (  FASMX  .LE.  FAS(I)  )  FASMX  «  FAS(I) 

IF  <  I  .LE  2  )  THEN 
FDRVT  -  DRVT 
OO  TO  640 
END  IF 

-  COMPUTES  FDRVT  FACTOR.  FREQUENCY  DEPENDENT  DURATON  OF  THE 

SOURCE  EXCITATION  IN  THE  RANDOM  VIBRATION  THEORY 

RATIO  «  FREQC I )  *  DRVT 

RATI03  -  RATIO  *  RATIO  *  RATIO 

RATI04  -  RATI03  /  (  RATIOS  ♦  1 .  0  /  3.  O  ) 

FDRVT  -  DRVT  ♦  RATI04  /  <  WI ( I >  *  DAMP  ) 

-  COMPUTES  OSCILLATOR  TRANSFER  FUNCTION.  H(FREQ. SFREQ) 

FREQ2  -  FREQ(I)  *  FREQ ( I ) 

FREQ4  -  FREQ2  •  FREQ2 

640  DO  700  J— 1 .  NSFR 

SFAS2  «  SFAS(J)  *  SFAS<J> 

SWI2  -  SWI(J)  •  SWI<J> 

-  SETS  H2  -  1.0  TO  COMPUTE  PEAK  ACCELERATION,  AMAX 

IF  (  I  .  EQ.  1  )  THEN 
H2  «  1.  0 
GO  TO  680 
END  IF 

-  SETS  H2  «  10/  SWI2  TO  COMPUTE  PEAK  VELOCITY,  VMAX 


I 


680 


700 


C 

C 

C 


C 

C 


IF  (  I  .  E0.  2  )  THEN 
H2  -  1.  0  /  8WX2 
00  TO  680 
END  IF 

FACT1  -  FRE02  -  8FREQ ( J )  •  8FREQ( J) 

FACT2  -  DAMP2  •  FREQ(I)  •  SFREQ(J) 

H2  »  FRE04  /  (  FACT1  #  FACT1  ♦  FACT2  *  FACT2  ) 

FUNXO(J)  -  8FAS2  *  H2 
FUNX2IJ)  -  FUNXO(J)  •  SWI2 
FUNX4( J)  -  FUNX2(J)  *  BWI2 

CONTINUE 

-  INTEGRATES  THE  OSCILLATOR  TRANSFER  FUNCTION  BY  SIMPSON'S  RULE 


CALL  SIHPS<  NSFR-1.  DSFR,  FUNXO.  SMO  ) 

SHO  -  SMO  *2.0 

ARVT  -  SORT (  SMO  /  FDRVT  ) 

CALL  SIMPS(  NSFR-1#  DSFR,  FUNX2.  SM2  ) 

SM2  -  SM2  *  2.  O 

PFZ  -  SORT (  SM2  /  SMO  )  /  PX2 

SNZ  »  2.  O  *  PFZ  *  DRVT 

IF  (  SNZ  .  LT.  2.  O  )  THEN 
FLAG  < I )  - 
SNZ  -20 
END  IF 

CALL  SIMPSt  NSFR-1.  DSFR.  FUNX4,  SM4  ) 

SM4  -  SM4  *20 

PFE  -  SORT (  SM4  /  SM2  )  /  PI2 

SNE  -20*  PFE  *  DRVT 

IF  (  SNE  .  LT.  2.  O  )  THEN 
FLAG(I)  -  '*' 

SNE  -20 
END  IF 

COMPUTES  OF ACT.  THE  RATIO  OF  AMAX  TO  ARVT. 


c 


$ 


& 


tv: 


$ 


720 


SION  -  SION  *  (-1.  O) 
CONTINUE 


C 

C 


OF  ACT  -  SR  HP  I  •  SUN 


ELSE 

OF AC  -  SORT (  2.  0  •  ALOO(SNZ)  ) 
OF ACT  -  OF AC  ♦  EULER  /  OFAC 
END  IF 


IF  (  I  .  EQ.  1  )  THEN 
PFEA  ■  PFE 
PFZA  -  PFZ 
SNEA  -  SNE 
8NZA  -  SNZ 
ARVTA  -  ARVT 
AMAX  -  ARVT  *  QFACT 
OO  TO  800 
END  IF 


IF  <  I  .  EQ.  2  >  THEN 
PFEV  -  PFE 
PFZV  -  PFZ 
SNEV  -  SNE 
SNZV  -  SNZ 

ARVTV  -  ARVT  *  981.  0 
VMAX  »  ARVTV  *  OFACT 
00  TO  800 
END  IF 


RSA(I)  -  ARVT  *  OFACT 
IF  <  RSAMX  .  LE.  RSA(I)  )  RSAMX  - 
RSV ( I >  -  RSA(I)  /  WI(I)  *  981  0 
IF  <  RSVMX  LE.  RSV(I)  >  RSVMX  » 


C 

C 

C 


800 


C 

C - 

C 


RSA ( I > 

RSV( I ) 

CONTINUE 

PRINT  1880.  IT-1.  AMAX.  ARVTA.  PFEA.  PFZA.  VMAX.  ARVTV,  PFEV,  PFZV 

IF  <  KEYS  .  EQ.  0  .  OR.  NIT1  EQ.  0  )  ©0  TO  860 
COMPUTES  SPECTRAL  RATIO  IN  THE  OSCILLATOR  FREQUENCY  REGIME 
DO  820  1*3, NFREQ 


IF  (  NPSV  .  ©T.  0  )  THEN 

SPRAT(I)  -  TSVI(I)  /  RSA ( I ) 
ELSE 

SPRAT ( I )  *  TSVI(I)  /  RSV ( I ) 

END  IF 


C 

C 


820 


IF  (  SPRMX  .  LE.  SPRAT(I)  )  SPRMX 
CONTINUE 


SPRAT ( I > 
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o  u  o  o  o  o  o  u 


c 


-  INTERPOLATES  SPECTRAL  RATIO,  SPRAT  AT  THE  INTEGRATION  FREQUENCY. 

SFREG(I) 

DO  840  I-l.NSFR 

CALL  INTERIM  FREQ.  SPRAT.  NFREQ.  SFREQ(I).  BSPRAT  > 

SPRATS ( I )  -  BSPRAT 
840  CONTINUE 

860  CONTINUE 


-  PRINTS  AND  WRITES  RVT  SPECTRAL  VALUES  OF  THE  LAST  ITERATION 

PRINT  1900 
PRINT  1200,  TITLE 
WRITE(8,  1100)  TITLE 
WRITE (8, 1100)  TITLE 
WRITE (8. 1920)  NFREQ— 2,  DAMP.  NIT1 
PRINT  1940,  DAMP*100.  O.  NIT1 
PRINT  1930 
C 

IF  (  NPSV  .  OT.  O  )  THEN 
PRINT  1960 
ELSE 

PRINT  1970 
END  IF 

C 

NLINE  -  10 
DO  870  1-3. NFREQ 

PRINT  1980.  1-2.  FREQ(I).  FAS(I).  RSA(I>.  RSVII),  TSVI(I), 

&  SPRAT  < I ) .  PER  < I ) .  FLAG ( I ) 

WRITE(8. 2000)  PER ( I ) ,  FAS(I>,  RSA(I),  RSV(I).  TSVI(I),  SPRATCI) 
NLINE  -  NLINE  ♦  1 
C 

IF  (  NLINE  .  CE.  60)  THEN 
PRINT  1200.  TITLE 
PRINT  1940.  DAMP* 100.  O,  NIT1 
PRINT  1930 
C 

IF  (  NPSV  .  OT.  O  )  THEN 
PRINT  1960 
ELSE 

PRINT  1970 
END  IF 
C 

NLINE  =  10 
END  IF 

C 

070  CONTINUE 
C 

CLOSE ( 8 ) 

PRINT  2020,  FASMX.  RSAMX.  RSVMX.  TSVMX,  3PRMX 
C 

IF  (  KEYS  EQ  0  )  00  TO  999 
C 

IF  (  NIT1  EQ  O  )  GO  TO  890 
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uu  o  u  u  u  o uo  o  o  <_>  o 


-  UPDATES  SFAS  WITH  THE  LAST  SPRATS 

DO  880  I-l.NSFR 
SFAS ( X  >  -  SFAS ( I )  *  SPRATS ( I ) 

880  CONTINUE 

-  COMPUTES  PHASE  FREQUENCY.  PFREQ(I)  AND  INTERPOLATES  SFAS  AT 

PFREQ(I)  FOR  INPUTTINO  TO  SUBROUTINE  TIME 

890  DO  900  I-l.NPFR 

PFREQ(I)  -  PFR1  +  DPFR  *  FLOAT (  I  -  I  > 

CALL  INTERPC  SFREQ.  SFAS.  NSFR.  PFREQ(I).  BFAS  > 

PFAS(I)  -  BFAS 
900  CONTINUE 

-  EXTRACTS  THE  PHASE  OF  THE  INPUT  ACCELERATION  TIME  HI8T0YR 

CALL  TIME(  -1.  0.0.  0.  0.0,  0,  CP.  SP,  OTH,  OTHMX  ) 

-  COMPUTES  OUTPUT  ACCELERATION  TIME  HISTORY  AND  RESPONSE  SPECTRUM 

OF  A  SINCLE-DEOREE-OF-FREEDOM  SYSTEM  <SDF>  FOR  NIT2  ITERATIONS 

DO  902  1=1 . NFREG 
PAA(I)  =0.0 
PRV(I)  =0.0 
902  CONTINUE 

C 

DO  912  IT=1 ,  NIT2+1 

C 

PAAMX  =00 
PRVMX  =0.0 
SPRMX  =00 
FACTOR  *  10 
C 

IF  (  IT  .  LT.  NIT2+1  )  THEN 

CALL  TIME<  O.  0.0.  O.  0.0.  O,  CP.  SP.  OTH.  OTHMX  ) 

ELSE 

CALL  TIM£<  O.  FC1.  NF1,  FC2,  NF2.  CP.  SP.  OTH,  OTHMX  ) 

C 

IF  (  FACTA  EQ  00  )  CO  TO  906 
C 

IF  (  FACTA  CT  0  0  )  THEN 
FACTOR  =  FACTA  /  OTHMX 
OTHMX  ■  FACTA 
ELSE 

FACTOR  *  AMAX  /  OTHMX 
OTHMX  *  AMAX 
END  IF 
END  IF 

C 

DO  904  1=1, NPT 

904  OTH(I)  =  OTH(I)  *  FACTOR 
C 

906  CALL  SPECT  <  NPT.  DT.  OTH,  PAA,  PRV  ) 


o  o  n  n  non  non  o  n  n  n  o  o 


C 


00  908  1-3,  NFREQ 

C 

FAS(I)  -  FAS!I>  *  8PRAT ( I ) 

C 

IF  (  PAAHX  .  LE.  PAA!I>  )  PAAMX  -  PAA<I> 

IF  !  PRVMX  .LE.  PRV! I >  )  PRVhX  -  PRV! I) 

C 

IF  (  NIT2  .EG.  0  )  GO  TO  908 
C 

IF  !  NP6V  .  OT.  0  )  THEN 

SPRAT!  I)  -  TSV I  <  I  )  /  PAAU) 

Pi  SE 

SPRAT! I)  -  TSV I ! I )  /  PRV!I> 

END  IF 
C 

IF  !  SPRHX  .LE.  SPRAT! I)  )  SPRHX  -  SPRAT! I) 

C 

908  CONTINUE 
C 

IF  (  IT  .  EQ.  NIT2+1  )  00  TO  912 

C 

DO  910  I-l.NPFR 

CALL  INTERB!  FREQ.  SPRAT.  NFREQ.  PFREQ!I>.  BSPRAT  > 

PFAS! I )  -  PFAS! I )  *  BSPRAT 
910  CONTINUE 
C 

912  CONTINUE 

-  UPDATES  FOURIER  AMPLITUDE,  FAS  IN  THE  OSCILLATOR  FREQUENCY  REGIME 

BY  BUTTERWORTH  FILTERING. 

FASMX  -0.0 

DO  914  1-3. NFREQ 

RE  -  CMPLX!  1.  O.  0.  0  ) 

HS  -  RE 

-  COMPUTES  THE  FIRST  BUTTERWORTH  FILTER. 

CALL  B WORTH!  FREQ! I).  FC 1 .  NF1.  HS  ) 

IF  !  FC1  .  LT.  0.  O  >  HS  -  1.0  /  HS 
RE  -  RE  *  HS 

-  COMPUTES  THE  SECOND  BUTTERWORTH  FILTER. 

CALL  B WORTH!  FREQ! I).  FC2.  NF2,  HS  ) 

IF  (  FC2  .  LT.  0.  0  )  HS  -  1.  0  /  HS 
RE  -  RE  *  HS 

FAS! I)  -  FAS! I)  *  CABS(RE) 

IF  !  FASMX  .  LE.  FAS! I)  >  FASMX  -  FAS! I) 

914  CONTINUE 

-  PRINTS  AND  WRITES  SDF  SPECTRAL  VALUES  OF  THE  LAST  ITERATION 
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c 


c 

PRINT  1200.  TITLE 
WRITEI9, 1100)  TITLE 
WRITE<9, 1100)  TITLE 
HR  I TE  ( 9.  1920)  NFREQ-2.  DAW.  NIT 
PRINT  2040,  DAMP*100.  0.  NIT 
PRINT  2090 
C 

IF  (  NPSV  .  OT.  0  >  THEN 
PRINT  I960 
ELSE 

PRINT  1970 
END  IF 

C 

NLINE  -  10 
DO  916  I-3,NFRE0 

PRINT  1980,  1-2,  FREQ(I),  FAS(I).  PAA(I>.  PRV!I>.  TSVI!I>. 

&  SPRAT! I).  PER! I) 

HRITE!9. 2000)  PER! I).  FAS! I).  PAA!I>,  PRV!I),  TSVI!I).  SPRAT! I) 
NLINE  »  NLINE  ♦  1 

C 

IF  !  NLINE  .  OE  60)  THEN 
PRINT  1200,  TITLE 
PRINT  20*0.  DAMP* 100. O.  NIT 
PRINT  2090 

C 

IF  !  NPSV  .  CT  O  )  THEN 
PRINT  1960 
ELSE 

PRINT  1970 
END  IF 

C 

NLINE  =  10 
END  IF 

C 

916  CONTINUE 

C 

CLOSE! 9) 

PRINT  2020.  FASMX .  PAAMX.  PRVMX,  TSVMX,  SPRMX 
C 

C  -  WRITES  OUTPUT  ACCELERATION  TIME  HISTORY 

C 

PRINT  1200.  TITLE 
PRINT  2200 

C 

WRITE! 1 1. 1 100)  TITLE 
WRITE! 11. 3000)  NIT 

PRINT  3100,  FI LEI 1 ,  NPT .  DT,  OTHMX 
WRITE! 11. 3200)  NPT.  DT,  OTHMX 
WRITE! 1 1. 3400)  !  OTH!I),  I-l.NPT  ) 

CLOSE! 1 1 > 

C 

IF  !  KEY6  .  EQ.  0  )  GO  TO  999 
C 

FACTA  «  FACTOR 
C 


! 


n  n  n  n  n  no  onoo 


C 


DO  918  I-l.NPFR 
PFAS(I)  -  PFASU)  *  981.0 
918  CONTINUE 
C 

00  TO  <  920.  930,  920  >  KEY 6 
-  COMPUTES  OUTPUT  VELOCITY  TIME  HISTORY 

920  CALL  TIM£<  1.  FC1.  NF1.  FC2,  NF2.  CP.  SP,  OTH.  OTHMX  > 

WRITE! 12,  llOO)  TITLE 
WRITE! 12, 4000)  NIT 

IF  <  FAC TV  .  EQ.  0.0  )  00  TO  924 

IF  (  FAC  TV  OT.  O.  0  )  THEN 
FACTOR  -  FAC TV  /  OTHMX 
OTHMX  -  FAC TV 
ELSE 

FACTOR  -  VMAX  /  OTHMX 
OTHMX  -  VMAX 
END  IF 

DO  922  I-l.NPT 

922  OTH(I)  -  OTH ( I )  •  FACTOR 

924  PRINT  4100.  FILE12.  NPT ,  DT.  OTHMX 
WRITE! 12. 4200)  NPT.  DT.  OTHMX 
WRITE! 12.  3400)  !  OTH!l).  I-l.NPT  ) 

CLOSE! 12) 

-  COMPUTES  OUTPUT  DISPLACEMENT  TIME  HISTORY 

930  CALL  TIME!  2,  FC1.  NF1 ,  FC2,  NF2,  CP,  SP,  OTH,  OTHMX  ) 

WRITE!  13,  HOO)  TITLE 
WRITE! 13. 3000)  NIT 

IF  !  FACTD  EQ.  0  0  )  CO  TO  934 

IF  !  r ACTD  CT.  0  0  )  THEN 
FACTOR  -  FACTD  /  OTHMX 
OTHMX  -  FACTD 
ELSE 

FACTOR  -  FACTA 
OTHMX  -  OTHMX  *  FACTOR 
END  IF 

DO  932  1*1, NPT 

932  OTH! I)  *  OTH! I)  *  FACTOR 
C 

934  PRINT  3100.  FILE13,  NPT,  DT.  OTHMX 
WRITE! 13. 3200)  NPT.  DT,  OTHMX 
WRITE! 13. 3400)  !  OTH!I),  1*1 . NPT  ) 

CLOSE! 13) 
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nonn 


999  PRINT  9990 


C 


STOP 


-  FORMAT  STATEMENTS 


1100  FORMAT  <  A80 ) 

1200  FORMAT <  '1 '»  //OX.  '***  ',ABO> 

1300  FORMAT ( //OX.  '***  COMMON  OUTPUT  FILE  NAME  :  '.  A80> 

1340  FORMAT! //OX.  '***  INPUT  PARAMETERS  ***', 

It  //9X«  '***  SOURCE  INFORMATION  ***', 

<i  //9X.  'STRESS  DROP .  ',F10.  3.  '  BARS'. 

«.  //OX.  'MEDIUM  DENSITY .  ',  F10.  3.  '  OM/CC  ' 

<>  //OX.  'EPICENTRAL  DISTANCE .  '.F10.  3.  '  KM', 

Si  //OX.  'SOURCE  DEPTH .  '.F10.  3.  '  KM', 

l<  //OX.  'SHEAR  WAVE  VELOCITY .  ',F10.  3.  '  KM/SEC 

«.  //OX.  'QUALITY  FACTOR,  Q .  ',F10.  3. 

t,  //OX.  'CONTROL  FREQUENCY  FOR  Q  FACTOR.  .  .  .  ',F10.  3.  '  HZ', 

«.  //OX.  'EXPONENTIAL  CONSTANT  FOR  Q .  '.F10.  3. 

«<  //OX.  'MOMENT  MAGNITUDE .  '.F10.  3. 

t,  //OX,  'SPECTRAL  DAMPINQ .  ',F10.  3,  '  X', 

Si  //OX,  'FMAX  CORNER  FREQUENCY .  ',  F10.  3.  '  HZ'. 

S.  //OX,  'FMAX  ORDER  NUMBER  .  ',  110. 

Si  //OX,  'EXPONENTIAL  FILTERING  CONSTANT.  .  .  .  ',F10.  3 

S.  //OX,  'KEY1  . '.110. 

&  //OX.  'KEY2 . 110. 

Si  //OX,  'KEY3 . ',110, 

S.  //OX,  'KEY4 . ',110, 

Si  //OX.  'KEY9 . ',110, 

S<  //OX,  'KEY6 . ',110) 


1400  FORMAT ( //OX.  '***  FUNCTIONAL  KEY  INFORMATION  ***  '  > 

1410  FORMAT ( //OX,  '***  KEY1  -  0  :  HALF-SPACE  RESPONSE  IS  COMPUTED') 
1420  FORMAT ( //OX,  '***  KEY1  -  1  :  SITE  RESPONSE  IS  COMPUTED') 

1430  FORMAT ( //20X,  '***  SITE  INFORMATION  ***', 

Si  //24X.  'DENSITY'.  6X.  'SV'.OX,  'DAMPING ',  OX.  'THICK',  OX,  'DEPTH 

!■  /24X.  '  (CM/CC) '.  3X,  '  (M/SEC) '.  OX,  '(X)',8X,  '<M)'.7X,  '(M)'. 

S.  /16X.  09< '—')»// 16X.  'SOIL  ',  9F10.  3, //16X.  'ROCK  '.2F10  3) 

1440  FORMAT < //OX,  '***  KEY2  -  O  :  NO  FREQUENCY  DEPENDENT  AMPLIFICATION 
S,  '  FACTORS  ARE  APPLIED') 

1460  FORMAT ( //OX,  '***  KEY2  -  1  :  FREQUENCY  DEPENDENT  AMPLIFICATION', 

Si  '  FACTORS  ARE  APPLIED  :  './/29X,  'FREQ.  '.6X,  'AMP. 

Si  /22X,  'NO.  '.  3X,  '  (HERTZ)  '.  4X.  'FACTOR  ',  /22X,  23 ('-'),  /> 

1480  FORMAT </20X.  10,  2F10  4) 

1000  FORMAT ( //OX,  '*#*  KEYS  -  1  :  ATTEN  -01/  SORT ( R )  -  ',F10  6, 

Si  '  FOR  R  >  100  KM.  ') 

1020  FORMAT ( //OX.  '***  KEY3  -  0  :  ATTEN  -  1.0  /  R  -  '.F10  6) 

1040  FORMAT ( //OX.  '***  KEY4  -  0  :  HERRMANN ' 'S  TERM,  HERR  -  0  O') 

1060  FORMAT ( //OX,  '***  KEY4  -  1  :  HERRMANN' 'S  TERM,  HERR  -  0.  00  *  R  - 

Si  F  1C.  4) 

1080  FORMAT ( //20X. '**#  SOURCE  EXCITATION  DURATION.  DRVT  -10/  FC'. 

Si  '  *  HERR  -  ',  F10.  4.  '  SECONDS') 

1600  FORMAT < //OX,  '***  KEYO  AND  KEY6  «  0  :  NO  SYNTHETIC  ACCELERATION', 
Si  '  TIME  HISTORY  IS  GENERATED') 

1620  FORMAT  <  //OX,  '***  KEYO  -  '.12.  '  :  SYNTHETIC  ACCEL.ERATION  TIME', 

S.  '  HISTORY  IS  GENERATED  BASED  ON  A  BUILT-IN', 


8.  '  ACCELERATION  TIKE  HISTORY ' > 

1640  FORMAT < //21X»  '***  INPUT  ACCELERATION  TIME  HISTORY  FILE  :  ABO. 

8.  //23X. 'FOR  MW  •  ',FS.  3.  '  AND  R  -  '.FlO.  3,  '  KM.  ', 

S<  //23X.  'WITH  NPIA  -  ',13.  '  AND  DT  -  ',F10.  3.  '  SEC.  ') 

1660  FORMAT ( / /3X»  '***  KEYS  -  '.12,  '  :  SYNTHETIC  ACCELERATION  TIME', 
l>  '  HISTORY  IS  OENERATED  BASED  ON  THE  INPUT', 

t>  '  ACCELERATION  TIME  HISTORY') 

1680  FORMAT ( /23X,  '*•#  ',A80) 

1700  FORMAT (21 3,  F10.  3,  A40) 

1710  FORMAT! /21X.  '***  WARNINO  :  ONLY  2048  POINTS  ARE  USED  IN  THE', 

<>  '  TIME  DOMAIN  COMPUTATION  !!!') 

1720  FORMAT < /21X,  '***  WARNINO  :  NP2  -  12  IS  USED  IN  THE  TIME  DOMAIN', 

I  '  COMPUTATION  INSTEAD  OF  '.12.'  !'.!') 

1730  FORMAT ( /21X. '***  FREQUENCY  RANOE  USED  IN  COMPUTINO  MODULUS  : ', 

6  //23X,  'NF  -  '.IS.'  .  FI  -  ',E13.  7.  '  HERTZ  AND  '. 

8.  'DF  -  ',  E13.  7.  '  HERTZ') 

1740  FORMAT ( //21X,  '***  NORMALIZING  FACTORS  FOR  OUTPUT  TIME  HISTORY  :  '. 

8c  //23X,  'FACTA  -  ',F10.  4,  '  0  ,  FAC  TV  -  '.FlO.  4.  '  CM/SEC'. 

8c  '  .  FACTD  -  '.FlO.  4.'  CM') 

1730  F0RMAT(//21X. '***  NO.  OF  ITERATIONS  -  ',12.'  :  WITH  NIT1  -  '. 

8t  12,  '  S<  NIT2  -  ',  12) 

1760  FORMAT < //21 X.  '***  DESIGN  (TAROET)  RESPONSE  SPECTRUM  :  ') 

1770  FORMAT  ( //29X,  'FREQ.  ',8X.  'SA'./22X.  'NO.  '.3X.  '  (HERTZ )'.  3X.  '(0"S)', 

8c  /22X.23(  '-'),  /) 

1780  F0RMAT(//29X.  'FREQ.  '.7X.  'SV'./22X.  'NO.  '.  3X,  '  (HERTZ  >',  2X.  '(CM/SEC)' 
Si  . /22X.  23( '-'). /) 

1790  FORMAT ( //29X. '***  OUTPUT  FILE  IS  :  ',A80> 

1800  FORMAT ( //3X,  '***  KEY6  -  ,  II.  '  :  OUTPUT  FILES  ARE  :  ') 

1820  FORMAT (  /20X,  1 1 ,  '  :  ',A80> 

1830  FORMAT ( //20X.  '***  FILTERING  PARAMETERS  FOR  OUTPUT  TIME  HISTORY  .  '. 


&  //23X.  'FC1  *  '.FlO.  4,'  HZ  .  NF1  -  ',13. 

8c  //23X.  'FC2  -  '.FlO.  4.  'HZ  ,  NF2-  '.I3> 

1840  FORMAT ( //3X, '***  COMPUTED  OUTPUT  PARAMETERS  ***', 

&  //3X,  'MEDIUM  DAMPING.  1/2Q .  ',F10.  3.  '  7.', 

S<  //3X.  'SOURCE  CORNER  FREQUENCY .  '.FlO.  3.  '  HZ', 

Sc  //3X,  'SOURCE  RADIUS .  ',F10.  3,  '  KM', 

V  //3X.  'SOURCE-TO-SITE  DISTANCE .  '.FlO.  3,  '  KM') 


1860  FORMAT ( / /3X,  '***  PEAK  VALUES  #**',//18X.  'PEAK  ACCELERATION'. 

S.  '  INFORMATION'.  14X.  'PEAK  VELOCITY  INFORMATION'. 

8c  /10X.  2(3X,  37(  '-')).  /6X.  'ITER.  '.  3X.  'AMAX'.6X.  'ARMS'.  7X. 

8.  'PFE  ',  7X,  'PFZ  '.  6X.  'VMAX'.6X.  'VRMS'.7X,  'PFE',7X. 

S.  'PFZ'.  /7X,  'NO.  ',  2(  3X,  '  (0  "S>  ' ).  2(6X.  '(HZ)').2X. 

8.  2(2X.  '(CM/SEC)  ').4X.  '(HZ)  '»  6X,  '(HZ)  '.  /6X.84(  '-'),/) 

1880  FORMATOX.  14.  IX.  2<2F10.  4.  2F10.  3)  ) 

1900  FORMAT ( ///3X,  '•  !  !  NOTE  .  *  -  FLAG  FOR  N  VALUE  <  2.  0  !  !  !  ') 

1920  FORMAT <13,  F3.  3,  3X,  'ITERATION  NO  '.12) 

1940  FORMAT ( /9X, '***  RVT  SPECTRAL  VALUES  AT  SPECTRAL  DAMPING  -  '.F3  2. 

8.  '  X  AND  AT  ITERATION  NO  ',  12) 

1930  FORMAT  </13X,  'FREQ.  '.  9X.  'FAS'.  12X,  'RSA  ' ,  12X.  'RSV',  12X.  'TSV'.  9X, 

S<  'SPECTRAL'.  3X,  'PERIOD' ) 

1960  FORMAT  ( 7X.  'NO.  ',  3X,  '(HZ)  '.  8X.  '(C.  SEC)  ',  10X,  '(G)  '.  9X.  '(CM/SEC)  '. 

V  10X.  '<©)',  11 X.  'RATIO'.  7X,  '(SEC)  ',  /7X,  98 (  '-').  /) 

1970  FORMAT  ( 7X,  'NO.  '.  3X,  '(HZ)  '.  8X.  '<C  SEC)  ',  10X.  '(G)  '.  9X,  '(CM/SEC)  '. 

Sc  7X,  '  (CM/SEC)  '.  9X,  'RATIO'.  7X,  '  (SEC)  /7X,  98(  >.  /> 

1980  FORMAT ( 7X. 13. FlO  4. 3E13  4, FlO  4, SX, A1 ) 

2000  FORMAT (6E 13.  7) 


2020  FORMAT (7X,  98<  '— ' }.  /7X»  'MAX.  VALUES  :  '.  3E13  4.  //) 

2040  FORMAT < /9X»  '**•  SDF  SPECTRAL  VALUE8  AT  SPECTRAL  DAMPINQ  -  '.F3 
L  '  X  AND  AT  ITERATION  NO.  12) 

2030  FORMAT  ( / 1 3X<  'FREO.  ',  9X,  'FAS  12X.  'PAA '.  12X.  'PRV'.  12X,  'TSV',9X, 
l  'SPECTRAL'.  3X,  'PERIOD') 

2200  FORMAT ( //3X«  '***  OUTPUT  TIME  HISTORY  INFORMATION  *** ' ) 

3000  FORMAT < 'OUTPUT  ACCELERATION  TIME  HISTORY  <0"8>  :  NIT  -  ',12) 
3100  FORMAT ( //9X,  '#**  '.A80,  //13X.  'NPT  -  ',13.  '  ,  DT  -  '.FlO.  4, 

*  '  SEC.  AND  AMAX  ■  ',  F10.  4,  '  O') 

3200  FORMAT!  13,  3X.  'DT  -  '.F10.  4.  '  SEC  .  AMAX  -  '.F10.  4.  '  O') 

3400  FORMAT (6E12.  3) 

4000  FORMAT ( 'OUTPUT  VELOCITY  TIME  HISTORY  (CM/SEC)  :  NIT  -  '.12) 
4100  FORMAT ( //9X,  '***  ',A80.  //13X,  'NPT  -  '.13.  '  .  DT  -  '.FlO.  4, 

«.  '  SEC.  AND  VMAX  -  '.  FlO.  4.  '  CM/SEC') 

4200  FORMAT  (IS,  3X,  'DT  «  ',F10.  4,  '  SEC  .  VMAX  -  ',F10.  4.  '  CM/SEC') 
3000  FORMAT < 'OUTPUT  DISPLACEMENT  TIME  HISTORY  (CM)  :  NIT  -  '.  12) 
3100  FORMAT  <  //9X,  '***  '.A80, //13X.  'NPT  -  '.IS,  '  .  DT  •  '.  FlO  4, 

S>  '  SEC.  AND  DMAX  -  ',F10.  4,  '  CM') 

3200  FORMAT  (IS,  3X<  'DT  -  '.  FlO.  4.  '  SEC  ,  DMAX  «  '.FlO.  4.  'CM') 

9990  FORMAT ( / / / 3X . '***  END  OF  PROORAM  ***'.////> 


FUNCTION  BINON (  L.  N  ) 


FUNCTION  BINON<  L. 


FUNCTION  BINON  CONFUTES  BINONI AL  COEFFICIENTS. 
BINON  -  N!/(L!  <N-L>  !  ). 

NOTE  :  L  MUST  BE  .  LE.  N 


BINON  -  1.  O 
DO  10  1-2.  N 

IF  <  I  . LE.  L  >  BINON  -  BINON  /  REAL ( I > 

IF  (  I  .LE.  <N— L>  )  BINON  -  BINON  /  REAL(I) 
BINON  -  BINON  •  REAL < I > 

10  CONTINUE 
RETURN 
END 


SUBROUTINE  BRUNEI  NFB.  NFE.  FREQ.  FA  ) 


SUBROUTINE  BRUNEI  NFB.  NFE.  FREQ.  FA  > 

SUBROUTINE  BRUNE  CONFUTES  FOURIER  AMPLITUDE.  FA  IO-SEC)  BASED  ON 
THE  METHOD  SUGGESTED  BY  J.  J*  BRUNE  (BS8A.  COL.  79.  SEPT.  1970). 
RELATED  PARAMETERS  ARE  CARRIED  OVER  FROM  THE  NAMED  COMMON  BLOCK. 
FADATA  FROM  THE  MAIN  PROGRAM. 


DIMENSION  FREQ!*)#  FA!*) 

COMMON  /FADATA/  KEY1.  KEY2.  C.  FC.  PIR.  Q.  FQ.  ALPHA.  FMAX.  N2. 
«.  CAP.  PICA  P 

COMMON  /PIDATA/  PI.  PI2 

COMMON  /AMDATA/  AFREG 1100).  AMFAC(IOO).  NAMP 
COMMON  /SMDATA/  SV1.  DAMP 1 .  THICK.  ARATIO,  RD.  RK.  CK 


AMPF  -10 
DO  300  I— NFB. NFE 
RATIO  -  FC  /  PREQ(l) 

-  COMPUTES  FREQUENCY  DEPENDENT  Q  FUNCTION 

IF  (  FQ  .  EQ.  0.0  OR.  ALPHA  .  EQ.  0.  O  >  THEN 
QFUNCT  «  Q 
ELSE 

QFUNCT  -  0  *  <  FREQ! I)  /  FQ  )  **  ALPHA 
END  IF 

ERATIO  -  PIR  *  FREQ(l)  /  QFUNCT 
FRATIO  ■  1.0  /  I  1.0  ♦  RATIO  *  RATIO  ) 

-  INTERPOLATES  FREQUENCY  DEPENDENT  AMPLIFICATION  FACTOR,  AMFAC  AT 

FREQUENCY,  FREQ(I) 

IF  (  KEY2  .  EQ.  1  )  CALL  INTERP ( AFREQ.  AMFAC.  NAMP.  FREQ ( I ) ,  AMPF) 

FA  <I)  «  C  *  EXPI— ERATIO)  *  FRATIO  *  AMPF 

IF  (  FMAX  .  EQ.  O.  O  .  OR.  N2  .  EQ.  0  )  GO  TO  100 

-  COMPUTES  BUTTERWORTH  FILTER 

FMFACT  -  SORT (  1.0  ♦  <  FREQ! I >  /  FMAX  )  **  N2  ) 

FA( I )  -  FA ( I )  /  FMFACT 

100  IF  I  CAP  EQ.  0.0)60  TO  200 

-  COMPUTES  NEAR-SITE  EXPONENTIAL  FILTER 

£  c 


FA<  I ) 


FA< I )  *  EXPI  -PICAP  *  FREQ 1 1 )  ) 


SUBROUTINE  BRUNE(  NFS.  NFE.  FREQ.  FA  > 


200  IF  (  KEY1  .  EQ.  0  )  00  TO  300 

-  COMPUTES  SITE  RESPONSE  AT  RECEIVER  DEPTH.  FAS  IN  THE  OSCILLATOR 

REOIME 

RK  -  PI2  •  FREQ(I)  /  8V1 
CK  -  RK  *  DAMP1 

CALL  SAMP (  ARDTH  ) 

FA< I )  -  FA(I)  *  ARDTH 

300  CONTINUE 


RETURN 


SUBROUTINE  B WORTH <  F»  FC,  NR,  HS  > 


SUBROUTINE  B WORTH (  F,  FC,  HP,  HS  > 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C  B WORTH  CALCULATES  THE  RESPONSE  OP  A  NP  POLE  BUTTERWORTH  FILTER 

C  UP  TO  AS  MANY  POLES  AS  THE  ARRAYS  8  AND  T  ARE  DIMENSIONED 

C 

C  F  -  FREQUENCY (HZ.  ) 

C  FC  -  THE  CORNER  FREQUENCY  OF  THE  FILTER 

C  NP  -  THE  NUMBER  OF  POLES.  NEQATIVE  FCR  HIGH  PASS 

C  HS  -  COMPLEX  RESPONSE  OF  THE  FILTER 
C 
C 

C  THE  FORMULA  USED  —  H(S)  -  1/<S-S1 >*<S-S2>  ...  (S-SK) 

C 

C  I*PI*< 1/2-M (2*K-1 >/<2*NP> ) > 

C  WHERE  SK  -  EXP  K  -  1,2.  ...  NP 

C 

C  S  -  I*(F/FC) 

C 

C  REF.  THEORY  AND  APPLICATION  OF  DIGITAL  SIGNAL  PROCESSING 

C  RABINE  t,  OOLD.  PAOE  227,  PRENTICE-HALL  1979 

C 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

COMPLEX  S ( 20 ) .  T ( 20 ) .  AS.  BK.  HS 
C 

N*IABS(NP> 

C 

C 

HS-CMPLX( 1.  . O  ) 

IF  (  N  .  EQ.  O  )  CO  TO  40 
IF  (  F  .  EQ.  O.  )  GO  TO  40 
SS-ABS(F/FC> 

DO  10  K-l.N 
AN  -  K 

AK  -  3.  141992694  *  <0.  9  +  (  (  (2.  *AN )  -1.  >  /  (2.  *N)  )  ) 

BK  ■  CMPLX (O. 0. AK) 

10  S(K)  *  CEXP(BK) 

C 

AS  *  CMPLX (0  0. SS) 

IF  (NP  LT.  0)  AS  *  1  /AS 
T ( 1 )  -  AS  -  S(l> 

IF  (N  .  EQ.  1 )  GO  TO  30 
DO  20  1*2.  N 

20  T ( I )  -  (AS  -  S( I > )  *  T(I-l) 

30  CONTINUE 
C 

HS  -  1  /  T(N) 

40  CONTINUE 

C 


RETURN 


SUBROUTINE  CMP MAX  <  KUO.  OT.  UQ.  PR.  w.  US.  W3.  WO.  D.  Z  ) 


SUBROUTINE  CMPMAX*  KUO.  DT.  UO.  PR.  U.  W2.  W3.  WO,  0.  Z  > 

C 

C  SUBROUTINE  CMPMAX  COMPUTES  RE8P0N8E  SPECTRUM  AT  LOW  FREQUENCY 

C 

C 

DIMENSION  U0<*>.  Z<*>,  C(3).  X<2.3) 

C 

00  10  1-1,3 
X<1.  I)-0. 

10  Z<I)-0. 

C 

FI  -  2.  *D/(W3*0T) 

F2  -  1.  /W2 
F 3  -  0*W 
F4  -  1.  /WO 
F5  -  F3*F4 
F6  -  2.  *F3 
E  -  EXP(— F3*DT> 

01  -  E*8 I N ( WD*DT ) 

02  -  E*COS <  WD*DT  > 

HI  -  WD*02  -  F3*01 
H2  »  WD*01  >  F3*02 
C 

DO  100  K  -  1.  KUO 
DUO  -  UO < K+ 1  >•  -  UQ(K) 

Z1  -  F2*DUC 

Z2  -  F2*U0<K> 

Z3  -  F1*DU0 

Z4  -  Zl/DT 

B  -  X(l.l>  +  Z2  — Z3 

A  -  F4*X (1.2)  ♦  F5*B  *  F4*Z4 
X(2,  1)  -  A *01  *  B*02  ♦  Z3  -  Z2  -  Z1 
X (2. 2)  -  A*H1  -  B*H2  -  Z4 
X <2.  3>—  F6*X<2.  2)-U2*X<2.  1 ) 

C 

DO  80  L-1.3 
C(L)-ABS<X(2.L> ) 

IF  <C(L)  . LT.  Z(L) )  00  TO  80 
Z(L)-C(L) 

80  X( 1. L)*X(2.  L) 

C 

lOO  CONTINUE 
C 

RETURN 

END 
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SUBROUTINE  COOL<  N.  X.  BIONI  > 


SUBROUTINE  COOL<  N.  X.  8I0NI  > 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

COOL  CALCULATES  EITHER  THE  FORWARD  OR  INVERSE  FINITE 
FOURIER  TRANSFORH  OF  A  COMPLEX  SERIES. 


NX 


F(J)  -  SUM  X(K)*EXP( I*8I0NI*2*PI*(K-1 >*( J— I >/NX> 
K-l 


J*l.  NX 


WHERE  NX  MUST  BE  AN  EXACT  POWER  OF  2. 


INPUT.  .  . 
N 
X 

SI  ON  I 


LOO (NX )  TO  THE  BASE  2 

COMPLEX  VECTOR  OF  DIMENSION  .  OE.  NX 

.  EO.  -1.  0  FOR  COMPUTATION  OF  FORWARD  TRANSFORM 

(FROM  TIME  DOMAIN  TO  FREQUENCY  DOMAIN) 
0  FOR  COMPUTATION  OF  INVERSE  TRANSFORM 
(FROM  FREQUENCY  DOMAIN  TO  TIME  DOMAIN) 


EQ.  4-1. 


OUTPUT. 

X 


THE  TRANSFORM  IS  STORED  IN  THE  POSITION  OF  THE 
ORIGINAL  SERIES.  THE  REAL  PART  CONTAINS  THE 
COSINE  SERIES  AND  IS  SYMMETRIC  ABOUT  THE  POINT 
2«*(N-1 )+l.  THE  IMAGINARY  PART  CONTAINS  THE  SINE 
SERIES  AND  IS  ASYMMETRIC  ABOUT  THE  POINT 
2»»(N-1)+1.  POINT  1  IS  FOR  ZERO  FREQUENCY  AND 
POINT  2**(N-1)+1  IS  FOR  THE  NYQUIST  FREQUENCY. 


NOTE.  .  . 

THE  COMPLEX  FACTOR 


1/NX 


1/2**N  IS  NOT  APPLIED. 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

COMPLEX  X(*>.  CARO.  CEXP.  CW.  CTEMP.  CPI 


CPI 

LX 

J 


CMPLX  <0.  0.3.  141396S338979*SI0NI > 
2**N 


-  1 


DO  30 
IF  (I 
CTEMP 


1-1. LX 
OT.  J) 
-  X(J) 
X( J)  ■  X ( I ) 
X(I>  -  CTEMP 
10  M  -  LX  /  2 


00  TO  10 


20 

IF  ( J  .  LE. 
J  ■  J  -  M 

M  -  M  /  2 

M) 

00 

TO 

30 

IF  (M  .  OE 

1  ) 

00 

TO 

20 

30 

J  -  J  ♦  M 
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SUBROUTINE  COGL<  N.  X,  SXONI  > 


40  I  STEP  ■  L  ♦  L 
DO  90  H-l.L 

CARO  -  CPI  *  <M— 1 )  /  L 
CW  -  CEXP(CARO) 

DO  90  I-M,  LX* I STEP 
CTEHP  -  CW  *  XCI+L) 

X(  I-H.)  -  XCI)  -  CTENP 
90  X(I>  -  X(I)  ♦  CTENP 

L  -  I8TEP 

IF  (L  .  LT.  LX>  00  TO  40 


RETURN 


SUBROUTINE  INTERS  (X.  Y»  N,  XX,  YY> 


SUBROUTINE  INTERS  < X. Y.  N. XX.  YY > 

SUBROUTINE  INTERB  LINEARLY  INTERPOLATES  ANY  P0INT8  BETWEEN  TWO 
SPECIFIED  POINTS.  IT  INTERPOLATES  YY  AT  XX  WITHIN  N  PAIRS  OF 
<X.  Y)  POINTS.  IF  XX  <  X(i>,  YY  18  8ET  TO  EQUAL  TO  1.0.  AND 
IF  XX  >  X(N>.  YY  IS  SET  TO  EQUAL  TO  1.0. 

DIMENSION  X(*>. Y<*> 

IF  (XX  .OT.  X <  1 ) )  00  TO  10 

YY  -  1.  O 

RETURN 

10  IF  (XX  .LT.  X (N) )  00  TO  20 
YY  -  1.  O 
RETURN 

20  DO  30  1-2,  N 

IF  (XX  .  LT.  X ( I ) )  00  TO  40 
30  CONTINUE 

40  YY  -  Y(I-l>«-(Y(I)-Y<I-l>>*(XX-X(I-i>>/(X(I>-X(I-l>> 


RETURN 


Ml 


VI 


SUBROUTINE  INTERP  ( X.  Y,  N*  XX#  VY) 


SUBROUTINE  INTERP  < X,  Y,  N,  XX.  YY) 

SUBROUTINE  INTERP  LINEARLY  INTERPOLATES  ANY  POINTS  BETWEEN  TWO 
SPECIFIED  POINTS.  IT  INTERPOLATES  YY  AT  XX  WITHIN  N  PAIR8  OF 
(X.  Y)  POINTS.  IF  XX  <  X(l).  YY  IS  SET  TO  EQUAL  TO  Y(l),  AND 
IF  XX  >  X < N) »  YY  IS  SET  TO  EQUAL  TO  Y<N) . 

DIMENSION  X<*>.  Y<*> 

IF  (XX  .  ©T.  Xd)>  QO  TO  10 
YY  -  Yd) 

RETURN 

10  IF  (XX  .LT.  X (N) )  00  TO  20 
YY  -  Y(N) 

RETURN 

20  DO  30  1-2.  N 

IF  (XX  .  LT.  Xd)  )  00  TO  40 
30  CONTINUE 

40  YY  -  Y(I-1)  +  (Y(I) — Y (I  — 1))*<  XX— X  d-1) )/<X(I)-X(I-l) ) 

RETURN 


on  n  r>  o  r> 


SUBROUTINE  SAMPt  ARDTH  ) 


SUBROUTINE  SAMP<  ARDTH  ) 

SUBROUTINE  SAMP  COMPUTES  SITE  AMPLIFICATION  AT  RECEIVER  DEPTH. 


COMMON  /SMDATA/  SV1,  DAMP  1 .  THICK,  ARATIO.  RD,  RK.  CK 


COMPLEX  CARD,  CATH,  ATH 
C 

RKRD  -  RK  •  RD 
CKRD  -  CK  *  RD 
C 

RKTH  -  RK  *  THICK 
CKTH  -  CK  *  THICK 

C 

CARD  -  CMPLX (  RKRD,  CKRD  ) 

CATH  -  CMPLX (  RKTH,  CKTH  ) 

C 

ARD  -  CABS <CCOS< CARD > ) 

ATH  ■  CCOS(CATH)  -  CMPLX<  0.0,  1.0  >  *  ARATIO  *  CSIN(CATH) 
ARDTH  -  ARD  /  CABS (ATH) 

C 

RETURN 

END 


64 


o  r>  o  o  o  o 


SUBROUTINE  SIMPS<  N.  DELTA.  FUNX.  AREA  > 


SUBROUTINE  SIMPS (  N.  DELTA.  FUNX.  AREA  > 

SUBROUTINE  SIMPS  INTEGRATES  INPUT  FUNCTION  BY  SIMPSON'S  RULE, 
NOTE  :  N  MUST  BE  AN  EVEN  NO.,  I>  FUNX(I>  MUST  HAVE  (N+l)  TERMS. 

DIMENSION  FUNX<*> 


SUM  -  FUNX  < 1 ) 

C 

DO  10  1-2,  N.  2 

10  SUM  -  SUM  ♦  4.  O  #  FUNX ( I )  ♦  2.  0  *  FUNX (1  +  1) 
SUM  -  SUM  -  FUNX ( N+l ) 

C 

AREA  -  DELTA  *  SUM  /  3.  0 
C 

RETURN 

END 
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SUBROUTINE  SPECK  NPT.  DT,  A,  PAA.  PRV  ) 


SUBROUTINE  SPECT(  NPT.  DT,  A,  PAA,  PRV  ) 

SUBROUTINE  SPECT  COMPUTES  A  RESPONSE  SPECTRUM  OF  A  SINGLE-DECREE- 
'FREEDOM  SYSTEM  IN  THE  TIME  DOMAIN  BY  SOLVINO  A  SET  OF  LINEAR 
SIMULTANEOUS  EQUATIONS  WITH  RECURSION  RELATIONSHIPS. 

DIMENSION  A<  * ) ,  PAA(*).  PRV<*>,  2<3) 

COMMON  /SPDATA/  DAMP,  NFREQ,  WK145),  PER < 143) 

KUO  -  NPT  -  1 

TTEST  •  10.  O  *  DT 

YY  -  SORT (  1.0-  DAMP  *  DAMP  ) 

DO  20  1*3.  NFREQ 
PR  -  PER ( I ) 

U1  -  Will) 

M2  «  W1  *  W1 
W3  -  HI  *  W2 
WD  -  W1  #  YY 

IF  (  PR  .  LT.  TTEST  )  THEN 

CALL  UCMPMX  (  KUO.  DT,  A.  PR,  Ul,  W2,  W3.  WD,  DAMP,  Z  ) 

ELSE 

CALL  CMPMAX (  KUO.  DT,  A,  PR,  Wl,  W2,  W3,  WD,  DAMP,  Z  ) 

END  IF 

PRV(I)  *  Wl  *  Z(l)  *  981.0 
PAA( I )  »  W2  *  Z( 1 ) 

20  CONTINUE 

RETURN 

END 


/-V-  /,  m 


o  o  r>  o  n  n  n  n  n  o  o  n  o  non  o  n  o  oo 


SUBROUTINE  TIME<  I FLA©,  FC1,  NF1 >  FC2,  NF2>  CP,  SP,  OTH,  OTHMX  > 


SUBROUTINE  TIME (  I FLAG,  FC1.  NF1.  FC2,  NF2,  CP,  SP,  OTH,  OTHMX  > 

C 

C  SUBROUTINE  TIME  COMPUTES  THE  FFT  OF  AN  INPUT  ACCELERATION  TIME 

C  HISTORY  AND  EXTRACTS  THE  PHASE.  IT  ALSO  COMBINES  AN  INPUT  MODULUS 

C  (BRUNE  FOURIER  AMPLITUDE)  WITH  THE  PHASE  AND  COMPUTES  AN  OUTPUT 

C  TIME  HISTORY  (ACCELERATION,  VELOCITY  OR  DISPLACEMENT). 

C 

C  THI8  SUBROUTINE  CALLS  A  SUBROUTINE.  COOL  (AN  FFT  ALOARITHM)  TO  DO 

C  THE  TRANSFORMATION,  AND  CALLS  A  SUBROUTINE.  B WORTH  (A  BUTTERWORTH 

C  FILTERING  ALOARITHM)  TO  DO  THE  FILTERING  FOR  COMPUTING.  OUTPUT 
C  TIME  HISTORIES. 

C 

C 

DIMENSION  CP(»).  SP ( * ) ,  OTH(«) 

C 

COMMON  /PIDATA/  PI,  PI2 
C 

COMMON  /TMDATA/  NPIA.  NPT ,  NPT2.  NP2,  NPFR,  DT. 

S>  A(4100).  PFREQ < 2030 ) ,  PFAS(2030> 

C 

COMPLEX  XX (4100) ,  RR,  HS 


IF  (  I FLAG  . OE.  0  )  GO  TO  300 

-  COMPUTES  THE  MEAN  OF  THE  INPUT  ACCELERATION  TIME  HISTORY 

SUM  -00 
DO  120  I-l.NPIA 
120  SUM  -  SUM  ♦  Ad) 

A ME AN  -  SUM  /  FLOAT <  NPIA  ) 

-  REMOVES  THE  MEAN  FROM  THE  INPUT  ACCELERATION  TIME  HISTORY 

DO  140  I-l.NPIA 
140  A<I)  -  Ad)  -  AMEAN 

-  ASSUMES  3  X  OF  NPIA  FROM  THE  INPUT  ACCELERATION  TIME  HISTORY  FOR 

COMPUTING  THE  COSINE  TAPER  AT  BOTH  ENDS  OF  THE  TIME  HISTORY 

NCT1  -  NINT (  0.03  *  FLOAT <  NPIA  )  ) 

NCT2  -  NPIA  -  NCT1  +  1 
CTFAC  -  PI  /  FLOAT (  NCT1  ) 

-  COMPUTES  COSINE  TAPER  AT  THE  FIRST  AND  LAST  NCT1  POINTS  FOR  THE 

INPUT  ACCELERATION  TIME  HISTORY 

DO  200  I-1.NCT1 


200 

Ad)  - 

Ad  )  • 

(  0.  3  - 

0. 

3 

*  COS< 

CTFAC 

*  FLOAT (I ) 

220 

DO  220 
Ad  )  - 

I-NCT2, 
Ad  )  • 

NPIA 

<03- 

0 

5 

*  COS< 

CTFAC 

*  FLOAT  < I ) 

-  ADDS  TRAILING  ZEROS  TO  THE  CONDITIONED  INPUT  ACCELERATION  TIME 

HISTORY 
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WVx, 


I AijLuittl  hwf> 


,  Yj&x  tX  A.  -2. 


SUBROUTINE  TXME<  XFLAO.  FC1.  NF1.  FC2.  NF2*  CP*  SP.  OTH.  OTHMX  > 


K> 


i 


DO  300  I-NPIA-d,  NPT 
300  A(I)  -  0.  O 
C 

C  -  ENTERS  FFT  CALCULATION 

C  -  TRANSFORMS  THE  INPUT  ACCELERATION  TIME  HISTORY  TO  THE  FREQUENCY 

C  DOMAIN 

C 

DO  400  I-l.NPT 

400  XX(I>  -  CMPLX (  Ad).  0.0  > 

C 

CALL  COOL<  NP2*  XX.  -1.0  ) 

C 

C  -  COMPUTES  THE  COSINE  AND  SINE  OF  PHASE 

C 

DO  420  I-l.NPFR 
RB  •  CABSC  XX (I )  )  *  DT 
CP (I)  -  REALt  XX (I)  )  /  RB 
SP(I>  -  A I  MAO  <  XXd)  )  /  RB 
420  CONTINUE 
C 

RETURN 

C 

C - INITIALIZES  A-  DUMMY  NORMAL I Z INO  COMPLEX  VARIABLE.  HS.  AND  RR 

C 

C  -  ENTERS  INVERSE  FFT  CALCULATION 

C 

500  RR  -  CMPLX (  10.  00  ) 

HS  -  RR 

C 

C - XXd)  IS  THE  COMPLEX  FOURIER  SPECTRUM  ARRAY 

C  —  :i(l)  IS  ZERO  FREQUENCY  ELEMENT  (  SET  TO  ZERO.  MEAN  -00) 

C 

XX  <  1  )  -  CMPLX  <  O.  O.  0  0) 

C 

C  -  ADDS  PHASE  <  CP.  SP  )  OF  THE  CONDITIONED  INPUT  ACCELERATION  TIME 

C  HISTORY  TO  SCALE  THE  BRUNE  MODULUS,  PFAS 

C 

DO  600  1-2. NPFR 

C 

XX(I)  -  PFAS<  I  - 1  )  *  CMPLX  (  CPd).  SP  d  )  > 

C 

C  -  COMPUTES  THE  FIRST  BUTTERWORTH  FILTER 

C 

CALL  BWORTHt  PFREQd-1).  FC1.  NF 1 .  HS  ) 

IF  <  FC1  LT  0  O  )  HS  -  I  0  /  HS 

RR  -  RR  •  HS 

C 

C  -  COMPUTES  THE  SECOND  BUTTERWORTH  FILTER 

C 

CALL  BWORTH<  PFREQd-1  >.  FC2.  NF2,  HS  ) 

IF  <  FC2  LT  0  0  )  HS  -  1  0  /  HS 

RR  -  RR  •  HS 

C 

XXd  )  -  XX  (  I  >  •  RR 
C 
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i 


I 

1 


SUBROUTINE  TIME<  I FLAG.  FC1.  NF1,  FC2.  NF2.  CP.  8P.  OTH.  OTHMX  ) 


RR  -  CMPLX  (  1.0.  0.0  > 

MS  -  RR 

PWI  -  PI2  •  PFREO(I-l) 

IF  (  IFLAO  .  OE.  1  )  HS  -  (  CMPLX (  0  0.  -PWI  >  >  **  I FLAG 

-  DIVIDES  THE  SCALED  SPECTRUM  BY  COMPLEX  OMEOA  **  IFALO  : 

IFLAO  -  O  :  COMPUTES  OUTPUT  ACCELERATION  TIME  HISTORY 

IFLAO  -  1  :  COMPUTES  OUTPUT  VELOCITY  TIME  HI8T0RY 

IFLAO  -  2  :  COMPUTES  OUTPUT  DISPLACEMENT  TIME  HISTORY 

XX<I>  -  XX<I)  /  <  HS  *  FLOAT < NPT)  *  DT  > 

-  NPT  *  DT  IS  THE  NORMALIZING  FACTOR  FOR  FFT  SUBROUTINE  COOL. 

INPUT  ACCELERATION  TIME  HISTORY  OF  DELTA  FUNCTION  OF  1/DT 
YIELDS  SPECTRAL  DENSITY  OF  UNITY. 

J  -  NPT 2  -  I 

XX(J)  -  CONJO(  XX(I)  ) 

HS  -  CMPLX (  1.  O.  O.  0  ) 

600  CONTINUE 

-  COMPUTES  THE  INVERSE  FFT 

XX(NPFR)  -  CMPLX (  REAL<  XX(NPFR)  ).  0.0) 

CALL  COOL<  NP2.  XX.  1.  O  ) 

OTHMX  »  OTHtl) 

DO  700  I-l.NPT 
OTH(I)  -  REAL (  XX(I)  ) 

IF  <  OTHMX  . LE  ABS ( OTH ( I ) )  )  OTHMX  -  ABS(OTH(I>) 

700  CONTINUE 

RETURN 

END 


nnnn 


SUBROUTINE  UCMPMX<  KUO.  DT.  UO.  PR.  U.  W2,  W3,  WD.  D.  Z  ) 


SUBROUTINE  UCMPMX<  KUO.  DT.  UO.  PR.  W.  U2.  W3.  WD.  D.  Z  ) 
SUBROUTINE  UCMPMX  COMPUTES  RESPONSE  SPECTRUM  AT  HIOH  FREQUENCY 


DIMENSION  UO(«).  Z<*>,  C(3>.  X<2.3> 
C 

DO  10  1-1.3 
X(l.  I)-0. 

10  Z(I>-0. 

C 

F2— 1 .  /W2 
F3— D*W 
F4-1.  /WD 
F3— F3*F4 
F6— 2.  *F3 
C 

DO  100  K-l.KUO 

NS— INT < 10.  *DT/PR-0.  01)  +  1 

DDT— DT /FLOAT ( NS ) 

FI -2.  *D/W3/DDT 
E-EXP <-F3*DDT) 

01— £*SIN(WD*DDT) 

02— E*COS ( WD*DDT ) 

H 1 -WD*C2— F3*Q 1 
H2— WD*0 1  -HF3*02 

DUO- ( UO ( K+ 1 ) -UO ( K ) ) /FLOAT ( NS ) 

O-UO(K) 

Zl— F2*DU© 

Z3— F1*DUC 
Z4-Z1/DDT 
C 

DO  lOO  IS— 1 . NS 
Z2— F2*0 

B— X< 1. 1>+Z2-Z3 

A— F4*X  <1.2)  -t-F3*B-*-F4*Z4 

X (2.  1 )— A*01+B*02+Z3— Z2-Z1 

X (2.  2) — A*H1— B*H2— Z4 

X ( 2.  3)—  F6*X<2.  2)-W2»X(2.  1 ) 

C 

DO  80  L— 1.  3 
C  <  L ) — ABS  <  X ( 2. L ) > 

IF  <C(L)  .  LT.  Z  < L )  )  00  TO  80 
Z<L>— C(L) 

80  X< 1. L)*X<2. L) 

C 

0-04- DUO 

c 

100  CONTINUE 
C 

RETURN 

END 
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Source  and  Propagation  Parameters 


Western  United  States  (WUS) 

Eastern  United  States  (EUS 

f  (cgs) 

2.7 . 

. 2.5 

(km/sec) 

3.2 . 

. 3.5 

*max 

15.0 . 

. 40.0 

k(kappa ) 

C.O . 

. 0.0 

Q(f) 

300.0* . . 

. 500.0  (f)0-65 

tCT  (bars) 

50.0 . 

. 100.0 

M0(cgs) 

log  M0  =  1.5  Mw  +  16.1 

log  M0  =  1.8  M.  +  14.1 

Z(f)  amplification 
factors 

See  Table  2 . 

. 1.0 

G(R )  ** 

R"1 

R"1 

V1 

V1 

fc 

tC  fi1!  8.44  M0 

ter £/8.44  M0 

•since  there  are  several  proposed  frequency  dependencies,  <pwas  left  constant  for 
these  predictions 


••geometri cal  attenuation 


m 


I 

$ 


i 


«'w» 


i 


:S 


*s 


i 


& 


1 


va  *"<1/5,  are  assumed  to  be  equal 
$  -  frequency 


0,R  refers  to  average  crustal  properties  and  near  receiver  properties  respectively. 


*T_  «r.  V*_ 


"A'AWi 


iiim 


Single  Layer  Site  Parameters 


(km/sec) 


/•  (cgs)  (X) 


(dampi ng) 


TABLE  4 


***  M4.  NF  DATA 


BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  M  <  4  3,  R 
OROVILLE  AFTERSHOCK,  CALIF.,  8/16/73.  ML  -  4  0.  R 
1320  PTS.  ,  DT  *  0.010  SEC  .  AMAX  -  0  068  0 


<  30  KM 
-  10  3  KM 


M4  NF  DATA 
CDMC  STA  P08S00E 


•••  M4  FF.  DATA 


BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  M  <  4  3,  R 
OROVILLE  AFTERSHOCK,  CALIF  ,  8/16/73,  ML  -  4  0,  R 
1320  PTS  .  DT  *=  0  010  SEC  ,  AMAX  «  0  068  C 


>  30  KM 
•  10  3  KM 


M4  FF  DATA 
CDMC  STA.  P08S00E 


***  M3  NF  DATA 

BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  4  3  <  M  <  3  3,  R  <  30  KM  M3  NF  DATA 
OROVILLE  AFTERSHOCK.  CALIF  ,  8/6/73,  ML  -  4  9.  R  -  8  6  KM  CDMG  STA  K06N33E 
1364  PTS  ,  DT  -  0  010  SEC  ,  AMAX  -  0  103  0 


***  M3  FF  DATA 


BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  4  3 

< 

M 

<  3  3,  R 

>  30 

KM 

M5 

FF 

DATA 

COYOTE  LAKE  EQ  , 

CALIF. , 

8/6/79, 

ML  -  3.  9, 

R 

m 

29  4  KM 

USGS 

STA 

CL10S37E 

2847  PTS  .  DT 

«  0  010 

sec. . 

AMAX  -  0 

106 

G 

***  M6  NF  DATA 

BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  3  3 

< 

fl 

<  6  3,  R 

<  30 

KM 

M6 

NF 

DATA 

COYOTE  LAKE  EQ  . 

CALIF  . 

8/6/79, 

ML  -  3  9, 

R 

■ 

18  4  KM 

uses 

STA 

CL02S40E 

2683  PTS  .  DT 

-  0  010 

SEC  . 

AMAX  »  0 

110 

G 

***  M6  FF  DATA 

BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  3  3 

< 

M 

<6  3.  R 

>  30 

KM 

M6 

FF 

DATA 

COYOTE  LAKE  EQ.  . 

CALIF  . 

8/6/79, 

ML  ■  3  9, 

R 

m 

29  4  KM 

USOS 

STA 

CL10S57E 

2847  PTS.  .  DT 

-  0  010 

SEC.  . 

AMAX  -  0 

106 

0 

***  M7.  NF.  DATA 

BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  M  >  6.  3,  R  <  30  KM.  :  M7  NF  DATA 
SAN  FERNANDO  EQ  ,  CALIF  ,  2/9/71,  ML  -  6.  4.  R  -  24  3  KM  :  USOS  STA  J144N69W 
2300  PTS  .  DT  -  0  010  SEC.  .  AMAX  -  0  283  0 


***  M7  FF  DATA 

BUILT-IN  ACCELERATION  TIME  HISTORY  FOR  M  >  6  3,  R  >  30  KM.  M7  NF  DATA 
IMPERIAL  VALLEY.  CALIF  ,  10/13/79,  ML  -  6  6,  R  >  40.0  KM  USOS  STA  IV23N43E 

2833  PTS  .  DT  -  0  010  SEC  ,  AMAX  -  0  110  0 


Earthquake  Magnitude 


WES  :  ATTENUATION 
MAGNITUDE  SCALING 


i 

(Log  Ap) 

( R  <  15  km) 


LEGEND 

-  While  rn  US  (  (I  —  10  km  ) 

.  Joyner  *  Boorc.  1982  (  B  -  10  km  ) 

-  Scad  it  Idr  I  is.  Bock .  1982  (  B  -  10  km  ) 

-  Donovan,  1973  (  B  -  10  km  ) 

| - 1  Sooklhs  *  Hanks.  1978  (  0  <  B  <  10  km  ) 

f- - 1  Joynor  it  Boorc.  1981  (  0  <  B  <  15  km  ) 

Figure  1.  WUS  Magnitude  scaling  at  close  distances  for  several 
attenuation  relations.  Each  relation  is  calculated  at  a  distance 
of  10  km  however  the  distance  definitions  differ.  Seed  and 
Schnabel  (1982)  and  Donovan  (1973)  employ  a  closest  distance, 
Joyner  and  Boore  (1982)  use  closet  distance  to  the  causative 
fault,  while  the  Brune  definition  is  for  hypocentral  range.  The 
data  ranges  (shown  by  vertical  bars)  are  for  hypocentral 
distances  of  less  than  15  km.  See  Table  1  for  WUS  parameters 
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Figure  2.  Peak  particle  velocity  at  a  distance  of  10  kn:  vs  moment 
magnitude  (Mw)  for  the  RVT  model  compared  to  the  prediction  of 

Joyner-Boore  (1982)  for  WHS. 
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Figure  5.  Peak  particle  velocity  at  a  distance  of  10  km  vs  magnitude 
(mO  for  several  EUS  attenuation  relations:  HBB,  Hasegaaw  et  al . 
(1981);  M,  McGuire  (1984);  A,  Atkinson  (1984).  Figure  taken  from 
Atkinson  (1984).  Open  symbols  are  from  the  RVT  model  with  EUS 

parameters  (Table  1 ) 
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PEAK  GROUND  ACCELERATION 

Figure  6.  EUS  peak  acceleration  distance  scaling  for  the  RVT 
model  (open  symbols)  for  mb  =  5  and  a  source  depth  of  4  km. 
Figure  taken  from  Atkinson  (1984)  with  data  as  indicated. 
Solid  line  is  Atkinson's  (1984)  Eastern  Canada  attenuation 

relation  for  mb  =  5 
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Figure  7.  Plot  of  significant  duration  (5%  to  95%  Arias 
intensity,  Dobry  et  al . ,  1978)  for  the  synthetic  time 
time  histories.  WUS  scaling  is  used  (Table  1)  and  the 
Brune  modulus  is  calculated  for  R  =  10  and  R  =  50  km. 
Based  upon  the  magnitude  and  distance  selection  criterion 
(see  Section  6.0)  appropriate  time  histories  are  automat¬ 
ically  selected  and  their  phases  combined  with  the  Brune 
modulus.  Also  shown  is  the  empirical  curve  from  Dobry 
et  al . ,  (1978)  with  the  t  2o  lines.  Interestingly 
enough,  the  inverse  corner  frequency  relations  (log  Tc 
=  0.5  M  -  2.49,  from  Table  1)  is  nearly  within  the 

4  2o  line 


Figure  Set  8.  Plots  of  synthesized  acceleration,  velocity,  and 
displacement  time  histories  for  WUS  and  EUS  tectonic  environments.  The 
event  has  a  magnitude  of  7  (Mw  for  WUS,  for  EUS)  at  a  hypocentral  range 
of  10  km.  The  phase  was  extracted  from  a  recording  of  the  1971  San 
Fernando  (M|_  =  6.4)  at  a  hypocentral  range  of  24.5  km.  Plots  of  response 
spectra  are  also  shown  following  the  displacement  time  histories. 

Also  shown  are  acceleration,  velocity,  and  displacement  time  histories  in 
addition  to  the  response  spectrum  for  the  WUS  event  recorded  at  a  depth  of 
20  m  within  a  40  m  thick  soil  site  (Table  3). 
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Figure  Set  9.  Plots  of  response  spectral  scaling  results. 


A)  Plot  or  response  spectra.  Open  circles  are  from  Boore  (1983)  and 
represent  regressions  on  WUS  data  for  a  moment  magnitude  6  event  at  a 
closest  distance  of  10  km.  The  upper  solid  line  is  initial  RVT 
response  spectrum  for  a  moment  magnitude  6  event  at  a  hypocentral 
range  of  10  km.  Dashed-dotted  line  is  the  specified  target  response 
spectrum  scaled  to  0.125  g.  The  dashed  line  is  after  two  iterations 
using  the  RVT  calculated  response  spectra  to  scale  the  Brune  spectra. 
The  solid  lower  line  represents  the  final  two  iterations  which  employ 
time  domain  calculations  of  the  response  spectra.  Time  histories  are 
calculated  by  adding  an  observed  phase  (from  the  time  history  library) 
to  the  scaled  Brune  modulus. 

B)  Plot  of  Fourier  spectral  density  after  two  iterations  (dotted  line) 
and  after  four  iterations  (solid  line).  The  first  two  iterations  have 
used  the  RVT  calculated  response  spectra  to  scale  the  Brune  spectra. 
The  final  two  iterations  use  time  domain  calculations  to  evaluated  the 
response  spectra. 

C)  Solid  line  shows  effects  of  band-pass  filters  on  response  spectra. 

The  filters  are  fifth  order  Butterworth  with  corners  at  4  seconds 
(high-pass)  and  23  Hz  (low-pass).  The  filtering  is  done  after  the 
last  iteration.  Remaining  curves  are  the  same  as  those  in  plot  (A). 

D)  Resultant  acceleration,  velocity,  and  displacement  time  histories. 

The  response  spectra  for  the  acceleration  time  history  is  shown  in 
plot  (C)  solid  line.  Unnormalized  peak  acceleration  is  0.116  g. 
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Figure  Set  10.  Plots  of  response  spectra  and  acceleration  time  history 
normalized  to  design  value  of  0.125  g. 

A)  solid  line  is  the  time  domain  calculated  response  spectra  of  the 
normalized  acceleration  time  history.  Remaining  curves  are  the  target 
response  spectra  (dashed-dotted)  and  the  second  iteration  RVT  response 
spectra  (dotted). 

B)  Acceleration  time  history  normalized  to  0.125  g. 
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Figure  Set  11.  Plots  of  results  from  a  different  scaling  technique  fitted 
to  the  same  target  response  spectra.  Target  and  final  response  spectra 
along  with  acceleration,  velocity,  and  displacement  time  histories  are 
shown. 
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Figure  Set  12.  Plots  of  basis  time  histories  which  are  used  to  supply 
phases  to  the  Brune  spectrums.  Magnitude  range  and  distance  range 
criterion  for  matching  to  the  Brune  spectra  are  given  in  Table  4. 
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COYOTE  LAKE,  MAIN  SHOCK,  1979 
ML  =5.9,  R  =  1 8.  7  KM. 
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